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INGENIEUR-ARCHIV 


Zur Theorie der rollenden Reibung 


Von H. Bufler 


1. Ubersicht. Wegen ihrer technischen Bedeutung haben sich schon mehrere Forscher mit 
-denjenigen Fragen auseinandergesetzt, die in mittelbarem oder unmittelbarem Zusammenhang mit 
dem Vorgang der rollenden Reibung stehen. Einen geschichtlichen Riickblick hierzu hat L. Féppl* 
-gegeben. Trotz zahlreicher experimenteller und theoretischer Untersuchungen gibt es jedoch noch 
verschiedene Probleme, fiir die bisher noch keine befriedigende oder noch gar keine Lésung gefunden 
wurde. Wahrend fiir gleiche elastische Konstanten der beiden in Beriihrung stehenden oder auf- 
einander abrollenden Walzen die wesentlichen Fragen geklart sind, stehen diese im Falle verschie- 
dener elastischer Konstanten noch offen. 
Im folgenden wird eine Ubersicht der mit der rollenden Reibung zusammenhangenden Pro- 
bleme gegeben, wobei eine Einteilung in drei Hauptfalle vorgenommen wird. 


Fall 1: Beide Walzen werden durch die Kraft P zusammengedriickt (Abb. 1a). 
A. Bei gleichen elastischen Konstanten der Walzen tritt langs der Beriithrungsflache keine Schub- 
spannung auf, und es gilt die gewéhnliche Hertzsche Druckverteilung. 
B. Bei verschiedenen elastischen Konstanten ist eine Schubspannung vorhanden, die ihrerseits 
die Druckverteilung beeinfluBt. Zu unterscheiden ist hier zwischen dem reinen Haftfall und dem 
-gemischten Fall, bei dem das Beriihrungsgebiet sich in eine Haft- und in zwei Gleitzonen gliedert. 


a b c 


Abb. 1. Beanspruchung der beiden Walzen. 
a) durch die Normalkraft P; b) durch die Normalkraft P und die Schubkraft Q (kein Rollen); c) durch die 
Normalkraft P und die Schubkraft Q, wobei beide Walzen aufeinander abrollen (Rad 2 sei das Treibrad). 


a) Haftfall: Dieser wurde von K. Desoyer® untersucht. Allerdings enthalten seine Ausgangs- 
-gleichungen Fehler, wie ich in Ziff. 2 zeigen werde, so dal er kein richtiges Resultat erzielt. 

b) Gemischter Fall: A. Schafer? behandelte dieses Problem unter Zugrundelegung der Hertz- 
schen Druckverteilung, was jedoch streng genommen nicht zulassig ist. 


Fall Il: Beide Walzen werden erst durch die Kraft P zusammengedriickt und dann mit der 


Schubkraft Q belastet, wobei kein Rollen stattfindet (Abb. 1b). 


1 L, Féppl, Die strenge Lésung fiir die rollende Reibung, Miinchen 1947. 

2 K. Desoyer, Osterr. Ing.-Arch. 11 (1957), S. 146. “ek ' y 

3 A. Schdfer, Die Reibung beim Walzendruck bei verschiedenen Elastizitatsmoduln. Diss. TH Miinchen 
(1948). 
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A. Bei gleichen elastischen Konstanten gilt die Hertzsche Druckverteilung, denn die Schubspan- 
nungen haben hier keine Riickwirkungen auf den iibertragenen Druck, wie erstmals H. Fromm} 
nachgewiesen hat. 

a) Haftfall: Dieser (nur bei unendlich groBem Reibungskoeffizienten mégliche) Fall ist relativ 
einfach zu behandeln (s. Gl. (55)). 

b) Gemischter Fall: Das Beriihrgebiet zerfallt hier in eine mittlere Haftzone und in eine vordere 
und eine hintere Gleitzone. 

c) Gleitfall: Die Lésungen fiir den Spannungszustand in den Walzen wurden unter Zugrunde- 
legung des Coulombschen Reibungsgesetzes angegeben von L. Féppl* und F’. Karas. 

B. Bei verschiedenen elastischen Konstanten unterscheidet man ebenfalls 

a) den Haftfall, 

b) den gemischten Fall und 

c) den Gleitfall. 


Fall III: Dieser Fall unterscheidet sich von Fall II dadurch, daB die beiden Walzen aufeinander 
abrollen (Abb. 1c). 

A. Bei gleichen elastischen Konstanten gilt wieder die Hertzsche Druckverteilung, die durch die 
Schubspannungen nicht beeinfluBt wird. 

a) Haftfall: Dieser ist nur bei unendlich groBem Reibungskoeffizienten méglich und wurde von 
L. Féppl* (S. 11) gelést. 

b) Gemischter Fall: Das Beriihrgebiet zerfallt hier in eine Haftzone auf der Anlaufseite und eine | 
anschlieBende Gleitzone. Die Lésung hat G. Heinrich® angegeben. 

c) Gleitfall: Hier gilt das unter Fall II, A, c) Gesagte. 

B. Fiir verschiedene elastische Konstanten liegen Untersuchungen von K. Desoyer® vor, und zwar | 

a) fiir den Haftfall, 

b) fiir den gemischten Fall und 

c) fiir den Gleitfall, der identisch ist mit Fall II, B, c). Allerdings sind die von Desoyer angege- | 
benen Ergebnisse wegen Fehler in den Ausgangsgleichungen nicht korrekt. 


Die vorliegende Arbeit beschaftigt sich in der Hauptsache mit den unter I, B, a) und II, A, a) | 
sowie III, B, c) genannten Problemen. Fiir die Druck- und Schubspannungsverteilung langs der _ 
Beriihrungsflache ergeben sich im Haftfall zwei gekoppelte singulare Integralgleichungen, die sich 
in zwei Punkten von den von K. Desoyer hergeleiteten Gleichungen unterscheiden und infolgedessen 
zu anderen Ergebnissen fiihren. Dariiberhinaus wird fiir die betrachteten Falle die tangentiale 
Randspannung innerhalb und auferhalb der Beriihrungszone ermittelt. 


2. Randdeformation. Da die Beriihrungslange 2 a sehr klein ist im Vergleich zu den Walzen- 
halbmessern, kann man beziiglich des Spannungszustandes zwei unendliche Halbscheiben zugrunde-- 
legen, die im Intervall —a S x S a mit der Druckspannung p und der Schubspannung q bean- 
sprucht werden. 


neyeleaelern 


_—— 


Abb. 2. Bezeichnungen an der Halbscheibe. 


Fir die Berithrbedingung und die Haftbedingung ist die Kenntnis der tangentialen 
Randdehnung (éx)y=o und der Ableitung der Verschiebungy nach x, dn/0x erforderlich. Die Ver- 
schiebung 7 zahle hierbei positiv ins Innere der Halbscheibe (wegen der Bezeichnungen, siehe Abb. 2) 


1 H, Fromm, Z. angew. Math. Mech. 7 (1927), S. 27 . 
2 L. Foppl, Forsch, Ing.-Wes. 7 (1936), S. 141. 
3 F’, Karas, Forsch. Ing.-Wes. 12 (1941), S. 266. 
4 Siehe FuBnote 1, Seite 137. 

> G. Heinrich, Osterr. Ing.-Arch. 4 (1950), S. 363. 

6 Siehe FuBnote 2 von S. 137. 
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Nach dem Hookeschen Gesetz gelten folgende Spannungs-Dehnungs-Beziehungen: 


é, = [o.—¥ (0, + a) (1) 
e, = 7 loy—¥ (0, +0,)1, (2) 
¢, = 7 [e,—(c. +0,)1- (3) 


Hierin bedeutet E den Elastizitatsmodul und » die Querdehnzahl. 
Setzt man einen ebenen Formanderungszustand voraus, so erhalt man wegen ¢, = 0 aus (1), (2), (3) 


1 
foe a) 91 +) oy], (4) 
1 
éy = [IL —*) 0, —v (1 +7) 0,]. (5) 
Fir einen ebenen Spannungszustand gilt dagegen wegen o, = 0 
1" 
Ex =F 16, —? oy] ? (6) 
1 
by = F ley — 62 : (7) 


Wie ein Vergleich der Ausdriicke (4), (5) mit (6), (7) zeigt, gehen die Formeln (4), (5) des ebenen 
Formanderungszustandes in diejenigen des ebenen Spannungszustandes (6), (7) formal dadurch 
uber, dai man »? = 0 setzt. Es geniigt daher, die Rechnung zunachst nur fiir den ebenen Form- 
anderungszustand durchzufiihren und erst die Ergebnisse auch auf den ebenen Spannungszu- 
stand auszudehnen}. 


Fir die tangentiale Randdehnung (e,), 9 der Halbscheibe erhalt man wegen (c,), 9 = — p 
aus (4) 

il 
(Ey =0 = 7 (CL —) (,)y=0 +L +¥) Pl: (8) 

Die Spannung (o,),—o betragt hierbei 

#1 
= = 2 £ gu) 
()y-0=—P@—= gb gg (9) 
=o 


(Das Zeichen ¢ besagt, daB das Integral wegen der Singularitat des Integranden bei u = x fiir x? < | 
als Cauchyscher Hauptwert zu verstehen ist.) Der Querstrich iiber einer Koordinate bedeutet hier 
wie im folgenden den mit der Strecke a (zunachst noch unbekannte halbe Lange des Beriihrgebietes) 
dimensionslos gemachten Ausdruck, also 

eee Fe reer 

Mes ies teers (10) 
Der erste Anteil in (9) rithrt vom Normaldruck p her (hydrostatischer Spannungszustand*); der 


zweite kommt von der Schubspannung q; denn fiir eine mit der positiven x-Achse gleichgerichtete, 
bei x = 0 angreifende Schubkraft Q (je Scheibendicke) gilt ® 


20 1 


(OZ) a ma x 


> 


woraus fiir eine kontinuierliche Lastverteilung durch Integration iiber den gesamten Lastbereich 


fiir den Aufpunkt x 
ail 
2 { qu) 7 
(o,)y <0 = —— ¢ re rae 
Si 


1 Streng genommen ist der sich tatsachlich in den Walzen einstellende Spannungszustand in der Umgebung 
der Beriihrungsflache dreidimensional. Um das Problem jedoch zweidimensional behandeln zu kénnen, werden 
als Vereinfachung die beiden Grenzfalle des ebenen Spannungszustandes und des ebenen Formanderungs- 

zustandes zugrundegelegt. Bei relativ diinnen Scheiben ist in Wirklichkeit der erstgenannte Fall zustandig, bei 
sehr breiten Walzen dagegen der zweite. 

2 L. Féppl, Drang und Zwang, Bd. 3, S. 37, Miinchen 1947. 

3 L, Foppl, a.a.O. 8. 29. 
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folgt. Gleichung (8) mit (9) ergibt somit 


a 
(c,)y-0 =— = | —» — 2”) p@) rea 6 wa 


—1 


(il 


Zwecks Ermittlung der Randverschicbung ny baw. der Ableitung @y/@x fiir die Halbscheib 
greifen wir auf (5) zuriick: 


e, = = — [(l—v) 0, (1 +») ,]. (12) 


Der Beitrag, der von p herriihrt, wird zweckmafig dadurch ermittelt, das man von einer im 
Koordinatenursprung angreifenden Normalkraft P (je Scheibendicke) ausgeht’, welche in der Halb- 


scheibe die Normalspannungen 
__2P By 
7 (x2 + y?)2 
und 
mde Mista! yi 
LG ts ma (x2 + y2)2 


zur Folge hat. 
Nach Einsetzen in (12) und Integration zwischen der oberen Grenze y = 0 und der untere 
Grenze y = h erhalt man die Randverschiebung y relativ zur Linie y = h 


x he 
Teall : 
a—Hn—=— +04) =| 


P 
(P 
We E 


x2 
Nimmt man him Vergleich zu x sehr groB, so folgt 


AOS —z3(¢ — 42) n(=) + 1 +?| 


Daraus bekommt man, wenn man auf eine kontinuierliche Lastverteilung iibergeht, 


ash 
= x —u\? ; — 
nf?) = ——* ¢ po la 52) n(® = *) +1+9|du. 
=) 
Der Differentialquotient 07{?)/dx ist somit 
4 
of?) 2 (1 —v*) pl) 4 
Ep tia yin Sapag pee (13% 
oo 
Um den Beitrag zu 7 zu ermitteln, der von der Schubspannung gq kommt, sei von einer im Be+ 
reich ¢ < x Sc konstanten, mit der positiven x-Achse gleichgerichteten Schubspannungsverteilung 
qo ausgegangen. Diese fiihrt zu den Normalspannungen? 
2 r q 
Cn aS = 3 (008 2 % — cos 2 a,) , 
und 
oy = — 3 (cos 2 a, — cos 2 a) , 
wobei 
Soa ea 2 y2 
ey = . 2 ol 2 — ss ye 
1 a Vy + (¢ eg ee cos 2 ay 1 y2+ (e+ x)?’ 
oe RNG =e cos 2.4 eee oe 
: P+ ee 


ist. Setzt man diese Ausdriicke in (12) ein und integriert iiber y, so bekommt man 


(a) = 0% |(] _y 992 Voce Mere | 
y aE | y YJ +2 +9) Ine eae 
wobei 


J =(ce —x) aretg —” — — (ce + x) arctg zie 


1 L, Féppl, a.a. 0. S. 61. 
2 I. Foppl, a. a. O. S. 39, 
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ist. Nimmt man als obere Grenze des Integrals y = 0, als untere Grenze y = o, so folgt 


RUE a 6 
nt) = 2 (L—v—22)a- c¢ fiir x>C, 
—c fir x<i—Ce 


und nach Differentiation 


en) 4 (L—v— 2) a for a! =. 6 
oe : (14) 
0 ing ei ek 


Bei einer veranderlichen Lastverteilung q = q(x) kann man eine Aufteilung der Belastung in 
einzelne Laststreifen q(x) dx = a q(x) dx vornehmen. Fiir jeden einzelnen Streifen gilt dann die 
Bezichung (14), die aussagt, da 07/dx proportional der Belastungsordinate an der betrachteten 


Stelle ist und unabhangig von den Belastungen auferhalb davon. Man kann daher allgemein 
schreiben 


on eet) 


oe Ez 4 q(x) . (15) 
Insgesamt gilt somit nach Uberlagerung von (13) und (15) 
a 
on a = 2(1— v? ee 
a= —E | —9— 20) aw +22” ¢ oy wu (16) 
241 


Die den Ausdriicken (10) und (16) entsprechenden Gleichungen fiir den ebenen Spannungs- 


zustand erhalt man, wie eingangs erwahnt, dadurch, daB man »? = 0 setzt: 


+10 
(e2y-0 = — 7 |(L—») p@) += Oa. (17) 
—1 
F j 2 fF 
zoe |-O—e@) +2 ¢ Pa (18) 
sl 


Gleichung (17) stimmt mit Gleichung (22) von K. Desoyer iiberein, 
Gleichung (18) weicht dagegen im Vorzeichen bei g von seiner aus 
Gleichung (3) mit (1), (2) resultierenden Formel ab, was auf einen dort 
vorhandenen Vorzeichenfehler zuriickgeht, der bereits in Gleichung (13) 
von G. Heinrich! enthalten ist. Ferner gewinnt Desoyer die fir den 


Ps ° 5 : Abb, 3. Zur Herleitung der 
ebenen Formanderungszustand maBgebenden Gleichungen aus denjenigen Hesthebediagunge 


des ebenen Spannungszustandes unrichtigerweise dadurch, da8 er statt E 

den Wert E/(1—v?) setzt. Die Unzulassigkeit dieses Schrittes zeigt 

der Vergleich der Formeln (17), (18) (ebener Spannungszustand) mit den Formeln (11), (16) (ebener 
Formanderungszustand). 


3. Beriihrbedingung. Zur Herleitung der Beriithrbedingung gehen wir von Abb. 3 aus.” Die 
elastische Anndherung der Kreismittelpunkte M, und M, ist, wenn 0 den (nachtraglich noch zu 
bestimmenden) Mittelpunkt des Beriihrgebiets bedeutet 


eels + + Ne + 


Ferner gilt aus geometrischen Griinden in erster Naherung 


AA =" Go + (7+ 2). 


im Ie 


Da fiir jeden Punkt des Beriihrgebiets 1 = konst ist, muB auch 0A/dx = 0 sein, woraus die Be- 
dingung 

é 1 = = Soe 

Ch +1) =—0(7 +7) Go +) fiir a? <1 (19) 


1 
ry To 


1 Siehe FuBnote 5 von Seite 138. 
2 Siehe L. Féppi, Drang und Zwang, Bd. 3, S. 67. 
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folgt. Entnehmen wir nun 07,/0x, und 072/0x, aus (16) (ebener Formanderungszustand) baw. (18) 
(ebener Spannungszustand) und beachten, dai gemaB Abb. 4 


x2 =X, Up =U, q2 => Pi — Pe 
OS re tS a UU ee 2 (i 18> 


ist, wobei der Index 1 sich auf die obere, der Index 2 sich auf die untere Walze bezieht, so erhalt 
man nach Einsetzen in (19) die Berithrbedingung in der Form 


amyl 


¢ PODS (EB ee fiir x2 <1 (20) | 


225 
Hierin bedeuten 


1 (il —==9F E, (1 —v?) beim ebenen Formanderungszustand, 
ee eit oa ee oul 
ry eS FE, + £, beim ebenen Spannungszustand. 
und 
1 2 2 

Ei 0 ee ere benen Formanderungszustand, 
CS aid E, (1 — ¥§) + E, (1 — vj) (22) 
2 


Elbe fe) mete eS 2) beim ebenen Spannungszustand. 
Bc 


Der fiir C mégliche Bereich ist 
TU Tt 
A= oe 


Falls Walze 1 (bzw. 2) als starr im Vergleich zu Walze 2 (bzw. 1) zu betrachten ist (E,/E, = 00 
baw. E,/E, = c0) und Walze 2 (bzw. 1) die Querdehnzahl null besitzt, so gilt C = 7/2 (baw. C = 


—z/2); bestehen dagegen beide Walzen aus demselben Material 
[ YT) (E, = E,, v1, = 2), dann wird C=0. 


Ly 


4, Haftbedingung. Wenn die beiden Walzen aufeinander 
abrollen, miissen in der Haftzone die sich beriihrenden Punkte 


v7) wery tt ‘rt tt ERAS dieselbe Geschwindigkeit haben. Diese Bedingung fiihrt, wie 
NALA 


von Desoyer gezeigt wird, in erster Naherung zu 


Ss S 
(a, = Exdy—0 9 HAW SEs) En) Oe (23) 
fiir x2< 1 


antreibt. Wesentlich dabei ist, daB der Schlupf S konstant 
Abb, 4. Randspannungen und Koordinaten an bleibt. 
Pesdon Saneiien (Walzen): Wir brauchen jetzt lediglich die Dehnungen aus (10) 
(ebener Formanderungszustand) bzw. (17) (ebener Spannungs- 
zustand) zu entnehmen und in (23) einzusetzen. Beachtet man dabei wieder, daB 


Ny =X, Uzg=—U, 2 — 4 > P2 == Ps 
Mp eee Us i) = = US Gh == Ole Par 


ist, so lautet die Haftbedingung 
41 


g di + Cp) =S" (im Haftgebiet) (24) 
—1 
wobei, falls Rad 2 im Sinne der Abb. 1c antreibt (diesen Fali wollen wir hier nur betrachten) 


g* — Ss ‘i ow is (1 — v3) + E, (1 —v?) beim ebenen Formanderungszustand, (25) 


Are E,+E#, beim ebenen Spannungszustand. 


C ist identisch mit dem in (22) angegebenen Ausdruck. 


| 
| 


je nachdem, ob Rad 1 oder Rad 2 (im Sinne der Abb. Ic) | 
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5. Der Gleitfall. Die beiden Walzen werden mit der Kraft P aufeinandergedriickt und das 
Antriebsmoment des Rades 2 sei so groB, da im ganzen Beriihrgebiet Gleiten eintritt (Falle II, B, 
c) und IIT, B, c))!. Unter Zugrundelegung des Coulombschen Reibungsgesetzes gilt dann 
| q(x) = + Mo p(*) (26) 
mit /o als Gleitreibungskoeffizienten. Das positive Vorzeichen in (26) ist durch die Antriebsrichtung 
_des Rades 2 gemaB Abb. 1c bedingt. Gleichung (26) ersetzt die Haftbedingung. 

Die Berithrbedingung (20) geht mit (26) iiber in eine singulare Integralgleichung 2. Art fiir p 
y Ns 
nC nop Ob ae fiir x2 
eet Po P(x) = a (%) + x) urx< 1. (27) 

1 
Die entsprechende Gleichung von Desoyer? weicht von Gleichung (27) durch das umgekehrte Vor- 
zeichen bei wf ab. Der Grund hierfiir ist der Vorzeichenfehler in seiner Formel fiir dy/0x. Durch 
einfache Umrechnung seiner Ergebnisse®, fiir die Druckverteilung p(x), die halbe Breite des Be- 
rihrgebietes a und die Lage des Mittelpunktes 0 (x9, siehe Abb. 3) ergeben sich mit den hier ange- 
wandten Bezeichnungen die korrekten Gleichungen 


p@) =" (eg| C49), (28) 


|/x2 ST we 1 +X 
py 
o=/5 FU— OB’ 2 
i=! =1— 20, (30) 
wobei 
B70 (0-270. 1) (31) 
0 
Die Schubspannung q(x) ist durch (26) mit (28) festgelegt: 
= Ba 1—x\? = 
= : 32 


AuBer der oben bestimmten Druck- und Schubspannungsverteilung interessiert der fiir die 
Beurteilung der Beanspruchung wichtige Verlauf der tangentialen Randspannung (6,),—o inner- 
halb und auferhalb der Beriihrungszone. Wegen (9) kénnen wir schreiben 


na GEM AR rs 2 q(u) a 3 
(o,)y=0 = — p(x) + = ¢ ee a fir xt ol (33) 
eal 
bzw. 
i 
ees q(u) es ( 
()y-0=+~ [48 fae: (3.4) 
et 


denn bei x2 > 1 ist wegen der freien Oberflache p = 0. Das obere Vorzeichen in (33) und (34) be- 
zieht sich hierbei auf Walze 1, das untere auf Walze 2. Mit (32) wird 
oa 


Ty 

UD) y= _ Ba Sf BV Le 

f erat = bo C3 ¢ tab 
1 


Das hier auftretende Integral laBt sich folgendermaBen umschreiben: 


+1 A 
¢ nase ee ea sefahmlarse Fy, 


aie ey) x—U 


1 Das ebene Problem eines starren Stempels beliebiger Gestalt, der auf dem Rand einer elastischen Halbscheibe 
gleitet, ist (nach ganz anderen Methoden) behandelt von N. J. Muskhelishvili, Some Basic Problems of the 
Mathematical Theory of Elasticity, Groningen 1953. 

2 Beim Vergleich der entsprechenden Formeln von Desoyer, hier (50), mit den hier abgeleiteten sind wegen 
der verschiedenen Festlegung der x-Achse die Indeces 1 und 2 zu vertauschen. 

3 K. Desoyer, a. a. O. Gleichungen (59), (58), (56). 
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Hierin bedeuten J, und J; die im Anhang abgeleiteten Integrale (4*) und (6*), falls man dort a = @ 
setzt.! Damit erhalt man aus (33) und (34) unter Beriicksichtigung von (28) und (31) 


—p@) (LF Hi C) + 7 Baw +129) qi ee de 
oe (35) 


(6,)y=0 = 2 “ © ao aes Pee 
+> Bajy|—@ +0 (EZ) +2 +1—26 fiir a> 1, 


wobei das obere Vorzeichen fiir Walze 1, das untere fiir Walze 2 gilt. Falls beide Walzen aus dem- | 
selben Material bestehen (C = 0, # = 1/2), so folgt aus (28) die bekannte Hertzsche Druckverteilung | 


p@) =—"/1—# (36) | 
und aus (35) 
yl—# operates he 


-l- 
2 2 
(o,)y-0 = + —Mo 4B isbolt (37) 
z——/e—1 fir @>1. 


Die Formeln (36) und (37) findet man auch sonst.” 

In Abb. 5 ist der Verlauf der mit der Bezugsspannung py = |/ P E (1/r, + 1/r2) dimensionslos 
gemachten Druckspannung p(x) fiir die in Tabelle 1 zusammengestellten Falle a bis d aufgetragen 
und zwar unter Zugrundelegung eines Reibungs- 
koeffizienten von (iy = 1. Die Kurve b gilt dann, 
wenn beide Walzen aus demselben Material be- 
stehen. Die Kurven a und c beziehen sich dagegen 
auf den Fall, daB eine der beiden Walzen starr 
ist (C= + 7/2). Zum Vergleich hierzu ist die 
Kurve d eingetragen; sie wurde fiir C = + 7/2 
bei Voraussetzung der Hertzschen Druckverteilung, 
welche sich fiir C + 0 nur im reibungsfreien Fall 
(Up = 0) einstellt, errechnet. Da die Konstante 
Cim Bereich — z/2 S CS a7/2 liegt und fy in der 
Praxis stets kleiner als 1 ist, stellen die Kurven a 
und ¢ praktisch kaum erreichbare Grenzkurven 
dar. In Wirklichkeit wird daher meist nur eine 
geringfiigige Abweichung von dem Spannungs- 


Abb. 5. Verlauf des Drucks p in der Beriihrungszone fiir die in ° 
Tabelle 1 angegebenen Fille a bis d bei vollstindigem Gleiten verlauf, der sich unter Zugrundelegung der Hertz- 
(Belastung nach Abb. lb oder le). schen Druckverteilung ergibt, vorhanden sein. 


Jedenfalls bleibt die GréBe des sich im allgemeinen 
etwas aus der Mitte verlagernden Druckmaximums in guter Naherung unverandert (vgl. die 
Extremkurven a und c mit d). 


Tabelle 1. Zusammenstellung einiger Spezialfélle. 


Fall E, E, Vy Vo Lo Cc o B 
A E | 00 |] | 1 es + 0,35 | —0,175 
b E E | | 1 0 0,50 0 

| 0 | 0 
c 0 E | 1 a 0,65 0,175 
d E (baw. 00) 00 (baw. E) | | Des > | 0,50 


1 Die mit einem Stern gekennzeichneten Gleichungsnummern beziehen sich auf den Anhang. 
* Siehe FuBnoten 2 und 3 von Seite 138. 
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Abb. 6 zeigt schlieBlich die tangentiale Randspannung (c,)y—o im Rad 2 fir die oben diskutierten 
Falle a, b und ec. 


(Oz)y=0 
% 


Abb. 6. Verlauf der tangentialen Randspannung (ox)y =0 im Rad 2 fiir die Falle a, b und c bei vollstindigem Gleiten 
(Belastung nach Abb. 1b oder Ic). 


6. Der Haftfall. Im Haftfall miissen die Beriihr- und die Haftbedingung entsprechend (20) 
und (24) erfiillt werden: 


sel ah] 
PO. di — Cz) = Ba +3), ¢ Sada + Cle) =S* furans -= | 
eo, Le 


Diese beiden gekoppelten Integralgleichungen fiir p und q lassen sich auf eine singulare Integral- 
-_gileichung 2. Art fiir die Funktion 


F(x) = p(x) + i q(x) (38) 


-zuriickfiihren: 
erent 
= [ LOGY se Si ] = = Fo ee 
FO eg p= ECB +A fie <1. (9) 


 Gleichung (39) ist vom Typ 
1 


De lan 
F@ ++ ¢ 2© di = Ce) farm? <1. (40) 
1 


Hi. Séhngen' gibt dafiir eine strenge Lésung in geschlossener Form an. 
Zunichst bekommt man durch Koeffizientenvergleich 


A=—it (41) 
mit |z/C| = 2 wegen (22) und 
- MIS) fer Baw, 4) Ba = 
C@) = 6 +6,% == im —iGe. (42) 


1 H, Sohngen, Math. Z. 60, (1954) S. 31. 
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Nach Satz 7 der Arbeit von Séhngen gilt : Ist G(x) in (—1, 1) eine Funktion der Klasse L? (p > 1) 
und werden solche Lésungen von (40) @ ain die in (—1, 1) einer Klasse L’(q > 1) angehéren, 
so gilt, wenn fiir J ++ +i 

1—iA 
1+iA 
ist: Ist A = iA, mit —1 <A, < 1, so hat die homogene Gleichung nur die triviale Lésung. Liegt 
A nicht in diesem Bereich, so wird die allgemeine Lésung der homogenen Gleichung gegeben durch 


ae 5 Eng e 1 Te x 
Fu) = — 0S =|» (44) 


= ei 2mm mit ee Ness) (43) 


15 
wobei 
+ 1 
C= Nelda, 
ey 
ist, und der Ausdruck 


ec 
Wx) G(x) A L 1—u/l—ul + x\* G(u) = : 
Ee) T+ 1+ a ¢ arte ee C1 urea, (45) 


liefert die allgemeine Lésung der inhomogenen Gleichung. Die Konstante 


fe 
C= f F(x) dx (46) 
= 
kann dabei beliebig vorgegeben werden. 


Da |x/C| = 2 ist, ist die Liésung (45) zustandig. Fiir « ergibt sich aus (43) und (41) 
wobei 
mt 
Gre 
Cem (48) 
ae | 
C 


ist. Die Konstante C errechnet sich nach (46 
ail 


C= [ F(x) de = — 
—1 —a 


Hierbei bedeuten P und Q die Normal- und Schubkraft je Walzendicke. 
Fiir die Funktion F(x) erhalten wir zunachst nach (45) mit Beriicksichtigung von (44) und 
(49) sowie (42): 


FG) =e — Tee we Tes) ae Oo +O — 9+ GLC + —D 2S} 


1—* 


P -QO\ snaa /l + x\% 
Slee ieee ee 


In (50) bedeuten J,, J, und J, die im Anhang angegebenen Integrale (4*), (5*) und (6*). Beachtet 


man die Beziehungen (47) und (48), so wird im einzelnen 
J, = (1—i2 8) yr — C2, 
1 = = = 
Jo= (5-28 +i2p))P—C, 


— 


und (38) zu 


e 


(p +iq) de =—(P +iQ). (49) 


—+ 


+ip 


| 
| 
t (51) 
h=—je—O@ +i cl, =3) (fiir #<)).| 


Die allgemeine Lisung fiir F(x) lautet dann wegen (50), (51), (41) und (42) 


B Ba Pred = 
(1 +46) 2S*h + > — Ba%z— Bax 


wey 


w +2 Baxp +2 Bap—S*)3], (52) 
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Vor der allgemeinen Diskussion dieser Gleichung sei noch auf den Spezialfall gleicher 


elastischer Konstanten der beiden Walzen eingegangen. Wegen C = 0 (d.h. fb = 0) geht 
(52) tiber in 


: 1 1 iP B = aa es 
FQ@) =p +iqg=+— {[e ae Baxyx— Baz) +i[T— svat. (53) 


Die Forderung, da® bei x = -L 


der Druck p verschwinden mu (Randbedingungen), verlangt: 
e PT aaa 
Xo = 0 und a = |/ 


woraus ; (54) 


als die Hertzsche Druckverteilung folgt. 

Die Schubspannung q sei zunachst fiir den Fall untersucht, da®B kein Rollen stattfindet, 
sondern bei vollkommener Haftung auf die mit der Kraft P bereits zasammengedriickten Walzen 
die Querkraft Q wirkt [Fall IT, A, a)]. Da unter diesen Umstanden keine Relativdehnung (¢,, — é,,) 
méglich ist, folgt wegen (23) S* = 0, und wir bekommen damit aus (53) 


a@) => = (55) 


y=0 


Diese Verteilung ware wegen der Singularitéten bei x = +1 nur bei einem unendlich groBen 
Reibungskoeffizienten méglich, so da sich in der Praxis neben einem mittleren Haftgebiet zwei 
seitliche Gleitzonen ausbilden. Falls die beiden Walzen jedoch aufeinander abrollen [Fall III, 
A, a)], ist ein Schlupf vorhanden, der im Haftfall ein reiner Formanderungsschlupf ist. Wegen 


(53) gilt hier 


Q 2 
Neigh lbalgs 2a 
OSS ji_z 


Fordert man, da8 auf der Anlaufseite (unter Zugrundelegung des Rades 2 als Antriebsrad ent- 
sprechend Abb. 1c), also bei x = 1, die Schubspannung null wird, so ergibt sich 


Ae aged 
und 
i ip x 
q(x) = — 2 ee (88) 


in Ubereinstimmung mit der von L. Féppl! angegebenen Formel. : 
Wir wollen uns nun dem Fall verschiedener elastischer Konstanten der beiden Walzen 
zuwenden und untersuchen, bei welchen Bedingungen hier Lésungen méglich sind. Aus physi- 


kalischen Griinden miissen wir fordern, daB bei x = + 1 der Druck null werden mu (Rand- 
bedingungen). Wegen (52) ist dies jedoch nur méglich, wenn 
Bate). (a) 
B 22 
> (1 +46’) —28*8 +—= Ba, (b) 
£42 Ba%mp = 0, (c) 
2Bap—S* =0.. | (d) 
Aus den Bedingungen (a) bis (d) erhalt ‘man 
; Ba 
%=0, Q=0, S*=2Baf8, P==—(1 +4). (57) 


Dieses Ergebnis besagt, daB bei vollkommenem Haften und unter Voraussetzung eines konstanten 
Schlupfes die Randbedingungen nur dann erfiillbar sind, falls keine Querkraft wirkt [Fall I, B, a)]. 
Die halbe Breite der Berithrungszone ergibt sich aus (57) zu 


Seer, 58 
o=Vazaaye: : Zs 


1 Siehe FuBnote 1 von Seite 137; (dort S. 11). 
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und fiir die Druck- und Schubspannungsverteilung erhalt man wegen (52) mit (57) die Formeln 
ae Ba SS 1+%x 
pa) = 7 1 —x* cos |8 In(; =| : (59) 


ws Ba see. Wh 35 
a) = | sin p In (i se =) (60) 
Die Ubertragung einer Normal- und Querkraft [Falle II, B, a) und III, B, a)] ist also bei voll- 
kommener Haftung fiir verschiedene elastische Konstanten der Walzen nicht méglich, selbst wenn 
man einen unendlich groBen Reibungskoeffizienten voraussetzt. Der Grund hierfiir liegt darin, 
daB die Ubertragung einer Querkraft auf unendlich groBe Druckspannungen an den Ecken fithren 
wiirde, was mit den Randbedingungen nicht vertraglich ist!. 

Wie aus (59) und (60) hervorgeht, wechseln p und q in unmittelbarer Nahe der Ecken x = + 1 
theoretisch unendlich oft das Vorzeichen. Der oszillierende Spannungsverlauf kann durch die 
wechselseitige Beeinflussung von p und q erklart werden und beruht auf der Voraussetzung eines 
konstanten Schlupfes. Da jedoch lings der Berithrungsflache kein Zug tbertragbar ist, folgt 
daraus, daB diese Voraussetzung nicht korrekt erfillt ist. Nach (59) und (60) liegt die erste Null- 
stelle von p im Falle C = z/2 bei 


x= ee fie 25-1 0ca1e 


also praktisch bei x = + 1, so da das sinusartige Verhalten auf die unmittelbare Umgebung 
der Eckpunkte, wo p und q fast schon auf den Wert null abgesunken sind, beschrankt ist. Das 
erste negative Maximum von p tritt fir C = 7/2 bei x = + [1 —3,7- 10-*] auf und betragt, wie eine 
einfache Rechnung zeigt, weniger als 1°{, des Drucks bei x = 0. Aus diesen Griinden kann man 
sagen, daB die Gleichungen (59) und (60) trotzdem brauchbare Lésungen fiir den Fall I, B, a) 
darstellen. 

Im folgenden soll der Verlauf der tangentialen Randspannung (¢,), —9 berechnet werden. Da- 
fiir gelten die Gleichungen (33) und (34). Wir gehen zweckmaBig von der aus (52) mit (57) folgenden 


komplexen Funktion 
ah Ba Aaa Le ee 


aus, welche fiir den betrachteten Fall zustandig ist, und are zunichst 


41 Nees 

RG) ea peas is 
¢ ne era ¢ —— du. (62) 
i) —l1 
Die Lésung fiir das vorstehende Integral findet man im Anhang, Gleichung (7*). Setzt man dort 
« =-—if und beachtet die Beziehung (48), so folgt 


= : C ——_= /1 — x\—i 3 
sbi reeggil bat BAG f fir @<1, 


J, = — yr? — C2(1 —i 2B) + wr— CP 


— 1\— ‘8 = 
ee y#—1 (| fir > 1, 
und man erhilt fiir den Imaginarteil von (62) 
1 Spel 
ga¢ Oa ¢ 1) oii 
x—uU x—uU 
af) i 
eae ag Vi—* —x2 cos B In(; =| fur = 
= Ba|2p — il (63) 
Ta —_ 1 sin E In E = || fiir x? Sale 
Nach Einsetzen von (63) in (33) und (34) unter Beriicksichtigung von (59) ergeben sich die Formeln 
4 = 2 
(6-0 = £5 Bap +p@) (FZ C—1) fir 2<1, 
(64) 


(y<0 = + 7 Bale — Tf sin pin(==))h fir 32> 1, 


wobei das obere Vorzeichen fiir Walze 1, das untere fiir Walze 2 zustandig ist. 


1 Vgl. hierzu H. Bufler, Z. angew, Math. Mech, 39 (1959) S, 218. 
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: Abb. 7. Druck p und Schubspannung q in der Beriihrungszone fiir die Fille a, b und c 
tei vollkommenem Haften und fiir den Fall d bei verschwindender Reibung (Belastung nach Abb. la). 


Abb. 8. Tangentiale Randspannung (Gx)y — im Rad 2 fir die Falle a,b und ¢ bei vollkommenem Haften 
Belastung nach Abb. 1a). 
g 


ta 


v8 


0 G5 1 
Abb. 9. Verhiltnisq/p in der Beriihrungszone bei vollkommenem Haften (Belastung nach Abb, la). 
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Abb. 7 zeigt den Verlauf der fiir die frither schon betrachteten Falle a, b und ¢ (siehe Tabelle 1; 
der Reibungskoeffizient 4) wird hier jedoch so groB angenommen, da® kein Gleiten méglich ist) 


berechneten Spannungen p und q, welche mit py = )/ PE (1/r, + 1/r,) dimensionslos gemacht sind. 
Zum Vergleich hierzu ist der Fall d eingetragen, der statt a und ¢ gilt, wenn keine Reibung vor- 
handen ist. Man sieht daraus, daB infolge der Reibung ein groBeres Druckmaximum auftritt als 
sich nach der Hertzschen Formel ergibt. Abb. 8 zeigt den Verlauf der tangentialen Randspannung 
(c,)y 0 im Rad 2 fiir die Falle a, b und c. In Abb. 9 ist schlieBlich das Verhaltnis p/q aufgetragen, 
und man erkennt daraus, daf der reine Haftfall nur bei einem unendlich groBen Reibungskoeffi- 
zienten méglich ware, so daf sich in Wirklichkeit neben einem mittleren Haftgebiet zwei seitliche 
Gleitzonen bilden. Auf dieses Problem, das eine gesonderte Untersuchung erfordert, soll jedoch 
an dieser Stelle nicht mehr eingegangen werden. 


7. Anhang: Zusammenstellung und Ableitung einiger Integralformeln. Aus einer Arbeit von 
H. Séhngen’ entnehmen wir folgende Formeln: ; 


+o 
iS We es2 — 1 1 ay ‘ Peer 
a | ate? = (— ta rd cee (0< ks<n, n ganzzahlig und positiv), (1*) 
a) 
1 “ ei—s)2 d wer - 
~¢ aug He gas (OMe see aul). (23) 


Hierbei kann s reell oder komplex sein; z und & sind reell. Ferner wird das anschlieBend abgeleitete 


Integral 
+00 
l > e(i—s)# d 25 Ome 0 (3%) 
Pal cata mrp coe oe) 


benétigt. Zum Beweis von (3*) betrachten wir voriibergehend z als komplexe Veranderliche. 
Der Integrand f(z) = h(z)/g(z) besitzt dann in der oberen Halbebene Pole erster Ordnung bei 
z,=F&+in(l+2n) mit n=0,1,2,..., falls 0< ¥is<1. Nach dem Residuensatz wird, 
wenn sich das Umlaufintegral titber die obere Halbebene erstreckt, 


; : h(z) ° oe sees : 
§ fe) de= 2b 3 Res dai 5 (MO) aad Sem aaaierver in ectinen A pat 


n=0 nO sings 


Durch Aufteilung des Integrationswegs erhalt man 


++ 00 mM 
Pfedz= [ f@dz+ | f(Reie) Ricivdg. 
—o 0 
reelle Achse R>o 


Halbkreisbogen 
Wie man leicht nachweisen kann, verschwindet jedoch der Halbkreisbogenanteil, womit (3*) 
bewiesen ist. 
Samtliche weitere bendtigte Integrale lassen sich auf die Grundformen (1*), (2*) ond (Gs) 
zurickfiihren. Mittels der Substitution 


= Go == |! 
US 
e + | 
wird 
ae! +a 
; we f eo? (L—«) 
an, = a dz, 
J i = aera | (@FIp 
und es folgt durch Anwendung von (1*) 
a0) 
ee he a 
i Fes Ott nee ol. ee (4*) 
il 
Ebenso erhalt man 
+1 4 ss ( +0 
a \e 2 —4) ee (et) 
——_—] udu = c 
i Fa) AOE ot ora ico | (e+ 1p 


* Siehe FuBnote 1 von Seite 145; (dort Gleichungen (19) fad (20)). 
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und nach Verwendung von (1*) 
1 


Ct Saah ee = 2 mt a? 
Ja= | (7G) He =— : (—1<Ra<)). (5*) 


—1 


ad i oe 
Beim Integral J; = ¢ van “/_ du hat man zwischen den Fallen x2 < 1 und x2> 1 zu unterscheiden'. 
=t 


Im Fall x? < 1 substituieren wir mit 
= ett = Cia 
a ag aig (a) 
und bekommen zunachst 

2 (1—a) 


+ co 
gs e 
J ———— (1 + e*) ¢ (f@+Dé —e) dz. 


Nach Partialbruchzerlegung des Integranden folgt 


+0 +o 
. et (L—4) q . ez 14) q b 
h= | Sayed See. (b) 
Im Fall x2 > 1 schreiben wir 
ee ee 
AE Cam ene (c) 
und man erhdlt in analoger Weise 
+0 + 00 
~ oz (1—a) * 2 (l1—a) 
ée€ e 
k= | Serpe | See (4) 


Wendet man auf die unter (b) und (d) stehenden Integrale die Formeln (1*), (2*) und (3*) 
an, so ergibt sich schlieBlich bei Beachtung von (a) bzw. (c) 


ee No ery = 
vee 1—(7F3) cose iin dag ean be 
Uaes ¢ ee ee Sees (6*) 
x—U sin 7 % x — toe Sp 
= 1-55 i) fur jaa 1 
mit der Einschrankung 0< Yia<l. 
Es soll nun noch das Integral 
pa 
nip be aro 
Tare ¢ ( + A] agree 
a4 
ausgewertet werden’. 
Im Fall x? < 1 wird mittels der Substitution (a) 
+0 
me g# B/2—2) ? 
Je ail) ¢ CAN TCT we 
und nach’ Partialbruchzerlegung des Integranden 
1 stj2—a) 1 getu2—a) 1 eti2—s) 
Cast me es Pi al aa 2 e MIS 
rite ial BONES Cok See = =. 
Jg= a soaps t2y fe | e+ 1 de LG e* — es : (e) 
Im Falle x2 > 1 erhalten wir mit (c) auf entsprechende Weise 
+a -+ co +0 
5 (1/2—a) Ey om Gare Se atta 0) 
Ae e on 2 e@& e ie" i ag ; (f) 
liane | ania 2 eae | Ere eect eet | ete 


1 Nur fiir x2 < 1 ist J, ein uneigentliches Integral. 
2 Nur fiir x2 < 1 ist J, ein uneigentliches Integral. 
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Die unter (e) und (f) stehenden Integrale lassen sich mit Hilfe von (1*), (2*) und (3*) berech- 
nen. Es folgt 


e& 


> [1+ sina eC? 49]. \ first 1, 


Dian if M Qn jire 
2 i errcene + sma) ot . 
ae [—1 + e—éC2+4)] fir <¢ 1. 
Beachtet man nun, da’ wegen (a) bzw. (c) 
fir x2< 1 
eo alae e211 — x? 
Loe 2 le 2 
wid furx? > 1 
é€ | 1+% e/2 x |/x—1 
le ee [ise = or eee ne 


so erhalt man: 


ltu 
im ————— /] — %\%", Pe =: 
1+3 +y1—F (ES) sin 7 a thie se << IL 
3 ; (7*) 
ce an Sabina af 
a eer B—1 (FF fir 2>1 


mit der Einschrankung — = <A Oe e ; 


8. Zusammenfassung. Nach einer Ubersicht der mit der rollenden Reibung zusammenhangenden 
Probleme werden fiir verschiedene elastische Konstanten der beiden Walzen die Berihr- und 
Haftbedingung sowohl fiir den ebenen Spannungszustand als auch den ebenen Formanderungs- 
zustand exakt abgeleitet. Die das Problem beschreibenden Integralgleichungen fiir den Druck p 
und die Schubspannung q innerhalb der Berithrungszone fiihren in den hier allein betrachteten 
Grenzfallen des vollkommenen Gleitens und vollkommenen Haftens auf geschlossene Lésungen. 
AuBer p und q wird die fiir die Beurteilung der Beanspruchung wichtige tangentiale Rand- 
spannung (0,)y—o0 berechnet. Die Ergebnisse werden diskutiert und fiir die méglichen Extrem- 
falle (eine der beiden Walzen ist starr; beide Walzen bestehen aus demselben Material) in Dia- 
grammen niedergelegt. Es zeigt sich, da8 im reinen Gleitfall gegeniiber der Hertzschen Verteilung 
im wesentlichen nur eine Verlagerung des Druckmaximums erfolgt, wahrend der reine Haftfall 
bei der Walzenpressung zu einer Erhéhung des Druckmaximums infolge der Reibung fiihrt. 


(Eingegangen am 21. Juli 1958) 
Anschrift des Verfassers: Dr.-Ing. Hans Bufler, Miinchen, Techn. Hochschule, Inst. fiir Techn. Mechanik 
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Ebenes Hochkantbiegen eines schmalen Balkens 
unter Beriicksichtigung der Verfestigung' 


Von Horst Lippmann 


1. Einleitung. Im Gegensatz zur Elastizitatstheorie sind die Grundgleichungen der Plastizitats- 
theorie bis heute keinesfalls unanfechtbar festgelegt. Selbst unter Voraussetzung eines plastischen 
Potentials (nach v. Mises*; vgl. z. B. auch Hill*) fiihrt die Wahl verschiedenartiger mathematisch 
zulassiger Formen desselben (wir erwahnen lediglich diejenigen von Tresca, Huber-Mises? und 
Coulomb, vgl. Shield*) auf z.T. ganzlich verschiedene Ergebnisse. Durchgefiihrte Versuche sind 
entweder zu einer endgiiltigen Entscheidung zu grob oder aber so fein, daB sie ohnehin aus dem 
Geltungsbereich der phanomenologischen Theorie herausfallen und mit dem physikalischen Modell 
der Kristallplastizitat (vgl. Kochendérfer®, Kréner®, Lippmann’) erklart werden miissen, welches 
seinerseits auf technische Formgebungsverfahren bisher nur in wenigen einfachen Fallen ange- 
wandt wurde (Lippmann’). Es gilt also, den auf technische Belange zugeschnittenen und daher 
stark vereinfachten Untersuchungen (wir verweisen beim Balkenbiegen z. B. auf die sehr interes- 
sante Arbeit von Hopkins®) theoretisch exakte Analysen bestimmter technischer Formgebungs- 
vorgange an die Seite zu stellen, die auch in Einzelheiten experimentelle Priifung zulassen. Wah- 
rend sich idealplastisches Blechbiegen geschlossen durchrechnen laBt?, muBte zur Beriicksichtigung 
der Verfestigung auf numerische Methoden zurickgegriffen werden (Proksa). Beim Hochkant- 
biegen fiihrt schon der idealplastische Fall auf keine geschlossenen Lésungen. Gaydon" bediente 
sich deshalb einer iterativen Naherung fiir nicht zu groBe Biegungen. Die vorliegende Unter- 
suchung erweitert Gaydons Ergebnisse auf verfestigenden Werkstoff und benutzt hierzu in Ab- 
schnitt 4 numerische Verfahren. Abschnitt 3 gibt auf Grund eines einfacheren plastischen Potentials 
allgemeine Bezichungen fiir die Geometrie des Vorganges — geschlossene vollstandige Integration 
ist dennoch unmdglich. Man erkennt vielmehr, auf welche verwickelten Ausdriicke selbst verein- 
fachte Ansatze fihren. Zu Beginn formulieren wir in Abschnitt 2 unsere Grundannahmen und 
Grundansatze. Naherungsbeziehungen fiir das elastische Werkstoffverhalten finden sich in Ab- 
schnitt 5, wahrend Abschnitt 6 einige grundsatzliche und tiberraschende Antworten zur Frage 
der Existenz und Eindeutigkeit der Lésungen erteilt: Bei gewissen Werkstoffverhaltnissen stellt 
sich eine Unvertraglichkeit im wesentlichen zwischen FlieBkriterium und den Spannungs-Gleich- 
gewichtsbeziehungen heraus, welche eine Unlésbarkeit der Grundgleichungen bedingt. Zwecks 
eingehender Analyse erweisen sich die z.T. etwas langatmigen Betrachtungen in 2b) und 6b) 


als unumgianglich. 


2. Voraussetzungen und Grundgleichungen. a) Grundannahmen. 

Voraussetzung |: Die plastische Formanderung ist groB gegeniiber der elastischen, so daB 
diese vorerst vernachlassigt wird. 

Demnach besteht die Alternative ,,starr-plastisch“. Linfliisse des elastischen Werkstoff- 
verhaltens beriicksichtigen wir naherungsweise in Abschnitt 5 (Riickfederung). 

Voraussetzung 2: Wahrend des ganzen Biegevorganges besitzt jeder Balkenquerschnitt Q 
eine Symmetriegerade g derart, da alle zu den verschiedenen Q gehdrigen g in einer Ebene £ 
liegen. 

Dies ist die exakte Formulierung fir ,;ebenes Biegen“. 


1 Dieser Aufsatz entstand am Lehrstuhl fiir Mechanik der Technischen Hochschule Hannover (Prof. 
Dr.-Ing. E. Pestel) mit Unterstiitzung der Deutschen Forschungsgemeinschaft. 

2 R.v. Mises, Z. angew. Math. Mech. 8 (1928), S. 161. 

3 R. Hill, The Mathematical Theory of Plasticity, Oxford 1950. 

4 R. T. Shield, J. Mech. Phys. Solids 4 (1955), S. 10. 

5 A, Kochendérfer, Plastische Eigenschaften von Kristallen und metallischen Werkstoffen, Berlin 1941. 

6 E. Kréner, Kontinuumstheorie der Versetzungen und Eigenspannungen. Erg. angew. Math., Hi. 5. 
Berlin/Géttingen/Heidelberg 1958. 

7 H, Lippmann, Acta met. 4 (1956), S. 298. 

8 H, Lippmann, Ing.-Archiv 26 (1958), S. 187. 

9 H.G. Hopkins, J. Mech. Phys. Solids 4 (1955), S. 38. 

10 F, Proksa, Diss. Hannover 1958, Stahlbau 28 (1959), S. 29. 

uF, A. Gaydon, J. Mech. Phys. Solids 1 (1952/53), S. 103. 
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Voraussetzung 3: Die Belastung des Balkens erfolge lediglich durch an den Endquerschnitten 


angreifende Momente. 

Insbes. sind alle anderen Begrenzungsflachen spannungsfrei. Tragheits- und Schwerekrafte 
werden vernachlassigt. 

Voraussetzung 4: Der Biegevorgang verlaufe in bezug auf Spannungen und Verzerrungen 
zentralsymmetrisch um eine gewisse (zu E senkrechte, raumfeste oder veranderliche) Achse mit 
dem DurchstoBpunkt P in E. 

Daher schneiden BalkenauBen- und Innenkante FE in konzentrischen Kreisbégen mit dem 
Mittelpunkt P (Abb. 1). Ferner miissen anfanglich ebene Querschnitte eben bleiben, da dies (aus | 
Symmetriegriinden) fiir den durch die Mittelebene PS herausgeschnittenen Querschnitt gilt. | 
Natiirlich zeigt sich erst im folgenden, ob das aufzustel- 
lende Grundgleichungssystem die Voraussetzungen 2 
und 4 gestattet. 


5 


VYoraussetzung 5: Die auf E wirkenden Span- 
nungen sowie ihre raumlichen Ableitungen verschwinden. 

In E herrscht ,,ebener Spannungszustand”. Er ist 
fir hinreichend ,,schmale“ Balken (deren Abmessungen 
senkrecht zu E gering im Vergleich zu den anderen 
Abmessungen sind) naherungsweise gewahrleistet. Dann 
bleibt aber auch die Querkriimmung der duBeren und 
inneren Begrenzungsflache ohne wesentlichen EinfluB: 


Voraussetzung 6: Die duBere und die innere Be- 
grenzungsfliche stellen Teile von auf E senkrechten 
Zylindern dar. 

Wir wahlen in E ein Polarkoordinatensystem mit dem 
(ohne Beschrankung der Allgemeinheit raumfest voraus- 
gesetzten) Punkt P als Pol und dem im Uhrzeigersinn positiv gezahlten Winkel 0 gegen 
die (ebenfalls raumfest vorausgesetzte) Richtung PS als Azimutalwinkel. Der Balken besitze 
die Anfangslange 1; seine Endquerschnitte schlieBen im betrachteten Augenblick den Winkel a 
ein (welcher den Biegezustand beschreibt): —oa = O <a. Alle auf den Anfangszustand « = Xo 
bezogenen GréBen mégen durch den Index 0 gekennzeichnet sein. Wegen Voraussetzung 4 und 
der Symmetrie um PS bilden die Koordinatenrichtungen Hauptspannungsrichtungen: o,, dg sind 
Hauptspannungen, €,, €g die zugehérigen Formanderungen. Partielle Ableitungen nach a werden 
durch Punkte bezeichnet. Ein Volumenelement besitze fiir « =a, die Koordinaten 0 = 0), 
r= Ry. Fira = o% liegt es bei O = Op a/ao, r = R(Ro, x). Es besitzt die Abmessung h = h(hg, x) 
senkrecht zu E und die ,,Radialgeschwindigkeit“ u(R,«) = R(R,, 4). Wegen des Werkstoff- 
zusammenhaltes fordern wir 


Abb. 1. Schnitt des Baikens mit E, 


Voraussetzung 7: u ist stetig in r sowie tiberall nach r differentiierbar. 

Das letzte gilt, damit é, [vgl. spater Gleichung (4b)] iiberall erklart ist. Es sei a(x) der Halb- 
messer der inneren, b(a) derjenige der auBeren Balkenkante und t(«) = b(a) — a(x) die Balken- 
dicke, sowie A(R, «) die am Volumenelement r = R bis zur Biegung « geleistete Arbeitsdichte. 
Soweit im folgenden (formal) benétigt, kennzeichnen wir die (nach Voraussetzung 5 verschwindende) 
auf E wirkende Normalspannung durch o,. Die zugehérige Verzerrungskomponente sei ¢,. 


b) Grundgleichungen. Die Misessche Theorie! des plastischen Potentials besagt die Existenz 
einer invarianten Funktion F(¢) des Spannungstensors o sowie einer skalaren Funktion A = A(r, a) 
mit 

F(o) =k (FlieBkriterium) , 


ened we F(o) (Stoffgesetz; i = r, O, h), (D 


wobei das Stoffgesetz tiberall dort gilt, wo die rechte Seite sinnvoll ist. Die sog. FlicBgrenze k = k(A) 
hangt stetig von der am betrachteten Volumenelement geleisteten Arbeitsdichte A(r,«) ab (vgl. 
Hill?, Abschnitt IT.3). 


Dene tubnete 2 auf Seite 153. 
2 Siehe FuBnote 3 auf Seite 153. 
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Wegen der Vernachlassigung von Volumenkraften und wegen Voraussetzung 5 lautet die einzige 
| Bedingung fiir das Spannungsgleichgewicht in E [siehe z. B. Szabé1, Gleichung (11.1)] 


0 
op (r 0,) — 06 . (2a) 


Aus den Randbedingungen o,(a) = o,(b) = 0 folgt 


i 09 dr == (0) 

| Setzen wir 
h/2 

Ko i o; dh (=F, @), 
l —h/2 
b 

so gibt { Kodr die auf Q wirkende Normalkraft, welche wegen Voraussetzung 3 verschwinden 
-mu8. Da die o; nur in E bekannt sind, liegt es nahe, K; = konst. o, zu setzen, wodurch wegen (2a) 


v 


b 
Kv eKdr 0 


-gewahrleistet wird. Allgemein kénnte man K;, fiir den Fall K 4 0 und Og = konst. ebenso aus- 
driicken. Es folgt 


b 
Ke 6, Olio, | hdr, 
also (auch im betrachteten Grenzfall K = 0) 


b 
fhdr 


K, (2b) 


oe 
: ba 


Oi 
; Kine zweite Méglichkeit besteht darin, (2a) fiir die K; niederzuschreiben (was exakt ist und 
_K = 0 gewahrleistet), dann aber fiir die 0; mittlere Spannungen K;/h zu benutzen. Dies bevorzugt 
2. B. Gaydon”; auch wir werden zum Vergleich mit seinen Ergebnissen i. a. auf diese Weise rechnen. 


‘Statt (2a) erhalten wir 


fa) : 
ap (r ha,) it Gian (3) 
Anders geschrieben, lauten (2a) bzw. (3) 
F 1 ; é 1 : 
na 109 Ol > (2’) zp (ht o,) = = [hog —ho,]. (3’) 


i 


Aus Griinden des Kraftegleichgewichtes in radialer Richtung verlangen wir tiberall 


Voraussetzung 8: Je nachdem, ob wir die linke oder die rechte Seite von (2’/3’) benutzen, 
sei o, baw. ho, stetig. 

Jedoch sind im folgenden Unstetigkeiten von og zugelassen. Daher mége (2’/3’) prazisiert 
werden: Uberall, wo die rechte Seite der Gleichungen stetig ist, mége die linke sinnvoll sein sowie 
(2’/3’) erfiillen und umgekehrt. 

In Anlehnung an Gaydon setzen wir naherungsweise 


(4a) 


Da wahrend der Verformung ebene Querschnitte in ebensolche tibergehen, ist die tangentielle 
Verschiebung bei Anderung von a um da von der Gestalt dv =r f(x, O) da. Z. B. bei Sommer- 
feld’ [Gleichung (4.28)], findet man 


ne, Szabo, Hohere technische Mechanik, Berlin/Gottingen/Heidelberg 1956. 
2 Siehe FuBnote 11 auf Seite 153. ; \ 
3 A. Sommerfeld, Vorlesungen iiber theoretische Physik II: Mechanik der deformierbaren Medien, Leip- 


zig 1949. 
11* 
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Da dies nach Voraussetzung 4 unabhangig von @ und dv = r da fiir O = « sein soll, ist f = O/a und 


__ ou 


: 1 A 
é=—+-, i, ==: (4b) | 


Wegen K = 0 muB gg langs r beiderlei Vorzeichen besitzen. Dann herrscht jedenfalls nahe bei | 

r =a Druck, nahe bei r=6 Zug: og <0 fir |r —al< oe; 09 >0 fir |r —b|<@ und hin- | 
reichend kleines 9. Wir setzen 
¢, = sup {7 |og = 0 und 0,0, fir aararye 


¢,==int {og = 0 fir 7 =r; 


gleich dem gréBten (kleinsten) Radius, unterhalb (oberhalb) dessen stets og << 0 und do S/o 
(vo > 0) gilt, und erkennen sofort a << ¢,<c, <b, ferner nach (2’), (3’) 6, << 0 fir acrSq | 
und c,<r<b. Das Trescasche bzw. Misessche FlieBkriterium lautet bei Unterdriickung der 
ohnehin verschwindenden Schubspannungsanteile 


F(o) = 0, — 99 =k bzw. : e | 
Gg = Vilar = ce F(o) = [(¢, — 46)? + (69 —9;)? + (6, —9,)7]? =k ; 
F(c) = 09 —0, =k bzw. wegen o, = 0 mit den Lésungen | 
09 =0, + V¥(.<rSb), tg =F —|p—to GarSc\e | (5) | 
wobei Y=k. Ser Se ats 
t= 5 +)/v—so (co<r<b), | 
| wobei positive Wurzeln gewahlt wurden und 
| Y=4h/y2 gilt. 


Die zugehérigen Stoffgesetze (bei Unterdriickung der verschwindenden Schubverzerrungen) sind 
&, =—te =A, é = 0) fir) asi cee. — 4 (0, ae ' 


ig = — é, =A; Ep 0 fiir Cr Ose Eo = tt (2 0g —G,) 5 (6) 


&, =" (—o,—9¢) mit p =A/k. 


Wir behaupten: Soweit iberall |o,|< Y (was bei kleinen a, also groBen a, b wegen (2’), (3’) 
stets zu erwarten ist), gilt c, = c, = c(x). 

Beweis: Die Werte c, und c, missen durch Unstetigkeiten von og gekennzeichnet sein. Die 
beiden Formen fiir og der rechten Seite von (5) gehen namlich héchstens fiir |o,| = 2 Y, /y3 >i 
stetig ineinander tiber. og = — Y der linken Seite kann nie, og = Y + o, nur fiir o, = — Y stetig 
das Vorzeichen wechseln. Wegen (4b) verlangt die Stetigkeit von u und g, in (6) 4 =0. Ebenso 
folgt A = 0; denn weil dort o, < 0 gilt, besitzt das Trescasche FlieBkriterium nahe c,, c, méglicher- 
weise die dritte denkbare Gestalt F(o) = o,—o, = —o, = Y, welche das Stoffgesetz é, = 0, 
— é, = & =A erzeugt. Bei c,, cy findet eine Vertauschung der drei méglichen Formen des Stoff- 
gesetzes statt, die wegen der Stetigkeit von ég stets A = 0 bewirkt. Wir finden jedenfalls mit 
(6), (4b) 

GF Ou(c;, &) 

a ae oa 
Es sei ¢, << cy, Dann besitzt (vgl. Abb. 2) u(r, x) wenigstens einen dazwischen gelegenen weiteren 
Schnittpunkt mit der Geraden —r/x. Offensichtlich gibt es einen von c, verschiedenen kleinsten 
solchen c;. Dort muB gelten 


4 ér  — Or 


u(cs, 4) =— =, RE A =| ai 
ee a 


Ce) 


Es folgt w # 0,440, & =0, also og =0,/2 oder |o,| = 2 Yy/V/3- auf den rechten Seiten von 
(5) und (6), wahrend die linken Seiten von (6) nicht gelten kénnen, so da8 ein FlieBkriterium 
lo,| = Y richtig sein miBte. Beides stellt einen Widerspruch dar, weshalb die Voraussetzung 
C, < €, falsch und der Beweis beendet ist. 

Wir beschranken uns zunachst auf kleine Biegungen derart, daB lo,|< Y. Hierfiir gibt uns 
(2’), (3’) in Verbindung mit einer der Seiten von (5) und (6) das vollstandige Gleichungssystem. 
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c) Transformation auf neue Koordinaten. Die bisher gewahlten Variablenr und a 
haben zwei wesentliche Nachteile: Erstens bewegt sich der Balken als ganzer im Polarkoordinaten- 
system (was den Vergleich verschiedener Phasen erschwert), zweitens werden. unsere Gleichungen 
fiir « = 0 (r = oo) singular. Dariiberhinaus bezieht sich « auf die willkiirliche und vom Problem 
unabhangige Ausgangslinge l. Umrechnung auf andere Ausgangslangen wird erschwert. 


ya rs 


Abb. 2. Méglicher Verlauf von u. 


Entgegen Gaydon transformieren wir bereits unsere Ausgangsgleichungen auf solche Koordi- 
naten, welche die oben genannten Nachteile vermeiden: 


sy = (Oy 1), 

b(a) (7) 
Pe (OS ea) 

b(x) — a(«) i 


Wir zeigen, da} y eine eigentlich monotone wachsende Funktion von a ist. Damit ist die Berech- 
tigung der ersten Substitution bewiesen. Die der zweiten sieht man sofort. 

Bekanntlich (vgl. Hill, Abschn. III.2.) gilt? 4» 20,4 = 0, also é SO fir aSrSc und 

é9 => Ofire <r<b. Ahnlich dem Beweis in b) schlieBen wir é < 0 fiir r = a, soweit |o,|< Y 

fir arc (bzw. é,>0 fiir r=6, falls |o,| << Y fir cS r<b): Andernfalls gilt namlich 

= 0 bzw. uw = 0 fiir r = a; das Gleiche nach b) fiir r = c. Wie dort (vgl. Abb. 2: c, = a, cg = ¢) 

erkennt man die Existenz eines Punktes r = d,a< d< c, wo ég = 0, jedoch wu # 0,4 ~ 0, was 
wie oben auf den Widerspruch |o,| = Y fiihrt. Nach allem ergibt sich also 


& ee ess, Cia iis 
é@ |, =b — €o eee 0 


u(b, «) u(a, &) 
a er mere se 


Wegen (7) ist dann 


= ; = 0. (8) 


oy _ a [u(b, a) a: u(a, @) 
Cea Sb a 


Damit ist die Berechtigung von (7) wenigstens im Falle |c,| < Y bewiesen. 
An Stelle von u(r, «) gebrauchen wir kiinftig die Funktion 


wey) = 5 |1 + ules) F |. 0) 


1 Siehe FuGnote 3 auf Seite 153. are: 
2 Im vorliegenden Fall 1a®t sich ohne Schwierigkeiten direkt beweisen: Fiir |c;| < Y erfiillt jedes Inte- 
gral unserer Grundgleichungen unter gegebenen Randbedingungen 4p = 0 bzw. A = 0. 
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wobei rechts r, « durch x, y ersetzt zu denken ist. An Stelle von R, R, benutzen wir X(Xo, y), Xo 
und gelangen zu folgenden transformierten, zunachst fiir y > 0 (a > 0) giiltigen Gleichungen: 

oe Ow. 

ery Ox” 

- y w—(U1—x) 

a aroe yee 


10 
__[ 1 eh} ey on 
Ep, —— Te a 0a ° 
2=~pL+(—y) w(l,y)— wy), | 
ferner 
OX w(X, y) — w(0, y) — X [w(1, vy) —w(0, y)] (11) 
ay yftd—y wvl.,y)— wy)” 
Da nach (6) é, + & + €, = 0 ist, so folgt 
422 py MO) Pa yao) = OO 12 
Bey yy et eee Om 
und schlieBlich 
oO ; _w(1, y) — w(0, y) = 
ay —") ya Fd —y) wy) — 00,90] (3) 


Dem Halbmesser r =c der neutralen Faser entspreche x = 2(y). Aus é9 =0 fiir x =z folgt 


wegen (10) 
w(2(y),y)) = 1 —2a(y). (14) 


Die Arbeitsdichte A errechnen wir aus 


Ow w — (1 — x) 
CASA pee Coney come Lae es ee a 
oy ay roy | 9 dy — yf[l1+(—y) wil, y)— w(, y)] * (15) 
0a 
Weiter ergibt sich 
pes iy, g bg 
Bera T=) ese) [% —o,], | Fy (h oD ear ae ay [h og —ho,] (16) 
(Gleichgewichtsbedingungen), 
0g = —Y fur) Oa Or 
® be uF, | de=>o—Y) ¥? 2 2 fur “OS 
Tea Os tures <a) | 4 
(Y=hk). | c= +/—t0 fir) 2 (1 
| = 
| (¥ = h/y2) 
(FlieBbedingungen nach (5)), 
[1d (hy eel yee a0. oh Fel em x) 0 | Wegen uw ~ 0 fiir x 4 z (Beweis 
( DD amet) oy a 1—y(1— x) >| wie oben mittels Abb. 2) 
Ow |dw  20,—oa w — (1 — x) 
— = < = Oyen as Nee 
a 9 (0S%x<2), Oy wd 6g = 1 = eee 
eee) 14+0—y) wily) wy) (8) 
én TY] mie | = 
Ue a) 
a Meche fyi Chea Ota ae 
yo (Gepa<ae lle (2 0@— 6, 1—y A — x) 
(CS). 


[Stoffgesetz nach (6) und (12)]. 
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a) Der ungebogene Stab. Fiir y = 0 bleibt die Koordinate x (im Gegensatz zu r) anschau- 
lich sinnvoll. Das Gleiche gilt fiir die w beschreibenden Differentialgleichungen unter (18), sofern 
wir deren Giiltigkeit im Falle y = 0 postulieren. Dies mége aus Stetigkeitsgriinden geschehen. 
Wegen der Symmetrie des geraden Stabes gilt 


= al! 
2(0) => iy 5 
~ nach (18) und (14) also in beiden Fallen 
| w(x, 0) = = : 


Postulieren wir ferner, y = 0 sei auch in bezug auf die Ableitungen von w und o nicht singular 
_ (dies wird sich im Falle des Trescaschen Kriteriums spater exakt zeigen), so folgt durch Differen- 
tiation der Differentialgleichungen fiir w von (18) nach y an der Stelle y = 0 


O2w 1 
om oe alle<*S}), 

aes l 

C2w 1 1 
ace a6 Aas ae es See 

be Mie ka) fir 5%, Ox dy |, — 9 eis | Ess )> 
Ww | 
wis eel Bey ate e  gee 

beso My fiir ES By ly 0 : A 


pel Sy Pn 1 
(m; = konst.). Damit dies fiir x = 5 
1 


stetig bleibt: m, — m, = a 


Fiir y > 0 streben Zahler und Nenner von (11) gegen Null. Anwendung der l’Hospitalschen Regel 
gibt als Grenzwert 


tee eas 2 =*(—x) 
Ox ma y¥=0 
ee 14 1 1 fir 0< x <1. ( (2) 
Dale 0 FU OiX—4) fr paxst. 


Beide Seiten sind (sogar fiir X = 1/2) stetig und stetig differentiierbar. Entsprechend erkennt man 
nach (13) 
ct 
oy 


t at | 


ee (a 13 
ects Bae (13a) 


3. Einige allgemeine Lésungsformeln. a) Formanderungsbeziehungen. Lediglich die zum 
Trescaschen Kriterium gehérigen linksseitigen Gleichungen gestatten teilweise geschlossene Inte- 
gration, aber auch diese nur in bezug auf x. Hierzu folgt aus (18) unter der Randbedingung (14) 


(a fue POS 4 = 2(4;) 
1 
oy) = 4 [Lz] y [1 —z)? + (L—x)?] (19) 
eR eee : z= far 2) i 
Ney (1 — 4) 
also wegen (13) 
| 62) = (lS 2 (20) 
oy Ey V2) Pe 12 y) G2)” 
wegen (11) 
eee) = fir 0<X<2(y), 
ax | 2it—yd—a) il + 12d—y)d—2)] (21) 
ye Cex ey ae) — 2] fir 2 X<1, 


20—ya— x) fi—yd— +L + 120—y d—a)] 
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Ppa DG 
ae — / fl—y@—xy f—yQ—404+120—y0—29ll 
j 0 fir xy ae 


b) Deutung und Vorschlage zur weiteren Integration. Spezialisierung von (20) auf 
y = 0 fiahrt genau auf (13a), von (21) auf (11a), was unsere dortige SchluBweise rechtfertigt. 

Die Gleichungen (19), (20), (21), (22) haéngen von der unbekannten Funktion 2(y) ab. Sie ge- 
wihren daher nur bei festem y einen Einblick in das Formanderungsverhalten. Wesentlich sind 
folgende allgemeinen Erkenntnisse nach (21) und (22): 

Der ganze ,,obere“ Balkenteil, welcher im Verlauf der Biegung noch nicht von 2(y) tiberstrichen 
wurde — wie wir unten aus Abb. 3 sehen, also mindestens die obere Balkenhalfte — besitzt seine 
urspriingliche Breite hy. Fiir alle anderen Balkenzonen findet Verbreiterung statt. Dagegen nimmt t 
eigentlich monoton ab. 

Ware z bekannt, so lieferte (21) eine explizite gewéhnliche Differentialgleichung erster Ordnung 
in X(X,, y) mit der Anfangsbedingung X(Xo, yo) = Xo. Sie ist teilweise sogar linear; im anderen) 
Teil diirfte sie kaum geschlossen integrierbar sein. Einfache (doch geschlossen kaum durchfiihrbare 
Quadratur liefert h nach (22). Dieses ganze Vorgehen wiirde durch die rechte Seite von (16) illuso- 
risch: Hier entstiinde eine Verkniipfung zwischen 2(y) und h. Statt dessen schlagen wir die linke 
Seite vor, welche in Verbindung mit (17) und vorgegebenem Y = Y(x) (insbesondere im idealplasti- 
schen Fall: Y = Y,) die Spannungen o,, dg sowie (als Schnittpunkt der von beiden Randern x = 0 
und x = | her errechneten Kurven o,) 2(y) zu ermitteln gestattet. 


c) Spannungen. Die Gleichungen (16) und (17) liefern fiir 0 < x S 2(y) 


fir” 0 ==, 


(22) 


a meats 
OV ee faery (Meee) J US ee 
andererseits fiir 2(y) Sx <1 
1 
Y 
Ong if l(a s des Gat 


Die Funktion z(y) geniigt der durch Gleichsetzen von (23a) und (23b) entstandenen Gleichung 


z 


f ¥@) de =[1—yQ—a] f age. (24) 


0 
Y = Y, = konst. fiihrt auf die transzendenten Lésungsgleichungen 
y+lny)=1—y, 
vy) =1—y [1 —2)]. 
Allgemein liefert der Mittelwertsatz der Integralrechnung 


1 1 
7 ¥(x) = acy dx 
Reerrerto by 


2 


(24a) 


f ¥(x) dx =Y’'s, 
6 


wobei Y, Y’ gewisse mittlere Werte von Y auf den Intervallen <0, z> baw. <z, 1) darstellen. Je 
kleiner y ist, desto symmetrischer verlauft Y um x = % = 1/2, desto genauer gilt also Y = Y’. 
Setzen wir diese Gleichheit annahernd fiir nicht zu groBe y voraus, so folgt wie oben (24a): z ist 
nahezu unabhiangig von Y; Y dient lediglich zur exakten F estlegung der Spannungswerte nach 
(22a), (22b) und (17) (betr. og). ; 


. Abb. 3 zeigt z(y) als Lésung von (24a) und demonstriert damit die Wanderung der neutralen 
aser, 
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d) Diskussion. Selbst die stark vereinfachte Theorie dieses Abschnittes laBt keine geschlossene 
Lésung zu. Sie fiihrt immerhin im Gegensatz zum allgemeinen Fall auf geschlossene gewéhnliche 
Differentialgleichungen. Doch verzichten wir auf deren numerische Integration, weil wir die auf 
dem Misesschen Kriterium beruhenden Gleichun- Z 
gen im nachsten Abschnitt numerisch auswerten. 
An dieser Stelle begniigen wir uns mit den oben 
genannten generellen Ergebnissen. 


4, Numerische Integration. a) Beschreibung 
des Verfahrens. Zur numerischen Integration 
legen wir von vornherein die Levy- Misesschen An- 
satze (rechte Spalte in Abschnitt 2) zugrunde und 
benutzen ein Differenzenverfahren, welches wie 


ie ied 
folgt arbeitet. a5 02 04 06 08 10 
Gegeben sei ein gewisser Zustand y mit den Abb. 3. 2(y) nech (24a). 
Koordinaten der Volumenelemente X(Xo, y) (wo- 
bei die Werte X, = X),, (« = 0, 1,..., N) 6, 


beim Anfangsstadium y = y, das Inter- % 
vall <0, 1) in N = 2 aquidistante Ab- 
schnitte teilen), den zugehérigen Breitenh 
und den Werten von Y. Ebenso sei die 
Balkendicke t bekannt. Die Schritt- 
weite Ay wird vorgegeben. Dann werden 
folgende Schritte zur Bestimmung des 
nachsten Biegezustandes durchgefiihrt: 

1. Erhéhung von y auf y + Ay. 

2. Berechnung der zu den X gehiri- 
gen Spannungeno,, dg mittels (16) und 
(17), wobei 0 durch A zu ersetzen und 
Ax = X(Xo i, y) — X(Xoj-1, ) von Stelle 
zu Stelle neu zu wahlen ist. Man 
findet den Schnittpunkt z beider Kur- 


venziige von o,. Genauigkeitskontrolle: 
1 


J (hoe) dx = 0. 
0 


3. Konstruktion von w und 


OES 


oe durch 
OX 


Integration der ersten Gleichung (18) 
mit der Anfangsbedingung (14). 
4, Errechnung der neuen Koordinaten 


X4+ Ss yaiictels (11) baw, (11a); 


5. Errechnung der neuen Breiten 


h+ = Ay nach der zweiten Glei- 
chung (18). 
Sea Ay mittels (13) baw. (13a). 


of ee + 5 Aynach (15)., 


8. Y mittels des vorzugebenden Ver- 


festigungsgesetzes Y = Y(A). 


b) Rechenprogramm und Ver- 
festigungsgesetz. Es wurde ein dem- 
entsprechendes Programm fiir den 
Magnettrommelrechner IBM 650 ent- 
worfen und erprobt, welches gestattet, 
wahlweise die rechte Seite von (5) durch 
die linke zu ersetzen. Dies stellt eine 0 025 05 0% 40 
Inkonsequenz gegeniiber der Theorie INN Hy. 


18 
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des plastischen Potentials dar, entspricht aber einem allgemeinen Brauch. Als Verfestigungsgesetz 


wahlen wir 
VY =Y. A\ PASS 
— (1 +r ( ) 


wobei y eine nicht negative Konstante ist und A die Werte 0, VRS Sr. : annehmen kann. Der Wertd 
hat also nichts mit der in Abschnitt 2 benutzten und in Abschnitt 6 wieder gebrauchten Funktion 4 
zu tun. 


bh 
15 Ro 
Yy=07 
125 | 
05 
93 
01 
10 ; 
Xe 
a ere a et 
479 02 04 06 08 70 
Abb. 4c. = 
12 
t 
t 
40 : 
y 
08 Easels 
0 Of 02 03 04 05 06 07 


Abb, 4d. 


In Zusammenarbeit mit dem Institut fiir Umformtechnik der Technischen Hochschule Hannover 
sind Versuche zur Priifung der Theorie geplant. Bei diesen wird Y(A) fiir die zu untersuchenden 


Werkstoffe experimentell (etwa mittels Zugversuches) festgelegt und an Stelle von (25) program- 
miert. 


c) Ergebnisse. Die folgenden Beispiele gehen simtlich von yy = 0, hy = konst., Y, = konst. 
aus und entsprechen einer Auswahl der von uns gerechneten Falle A = 1, 2,3; OS y»<5. Der 
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Wert » = 0 bedeutet idealplastischen Werkstoff (Abb. 4) und gestattet einen Vergleich mit den 
Naherungen Gaydons. In der Tat stimmen diese fiir kleine y fast vollstandig iiberein. Ledighich 
fiir groBe y sinken entgegen Gaydon 
sowohl Stabdicke t wie Biegemoment 


1 
Vit [vow dt (26) 
0 


ab [Betreffend x; vgl. 2 b)}. Das steilere 
Absinken von M bei Benutzung von 
Trescas Kriterium entspricht wiederum 
Gaydons Ergebnissen. 


Interessanterweise ist der Verlauf 
von fh praktisch fiir alle Rechnungen — 
unabhangig vom Grad der Verfestigung— 
derselbe. Ahnliches finden wir bei t. 
Die Spannungen variieren in ihrer GréBe, 
weniger im charakteristischen Verlauf 
(Abb. 5, 6, 7). Dabei bleiben die nach 
Tresca errechneten Spannungen und 
Momente wie bei Gaydon bis zu 20° 
unter denen nach Mises. 


Abb. 8 gibt einige Beispiele von 
Momentenverlaufen, Abb. 9 von Y-Ver- 


teilungen (nur fiir den Biegeparameter 0 EG 05 D5 70 
y = 0,7). Abb, 5a. 

d) Analyse und Vergleich der 
Ergebnisse. Die aufgezwungene Form- ty 


anderung y scheint die verdnderlichen 
Abmessungen hf, t unabhangig vom 
Verfestigungsgrad zu bestimmen, so da 
diese im wesentlichen durch Abb. 4 ge- 
gebensind. Ein experimenteller Vergleich 
laBt daher lediglich auf die Richtigkeit 
der Theorie, keineswegs auf die des 
Verfestigungsgesetzes schlieBen. Hieran 
andert die Benutzung des Trescaschen 
Kriteriums nichts. Dies gilt nach den 
Uberlegungen von 3c) (z unabhangig 
vom Verfestigungsgesetz) zwar auch bei 
konsequenter Benutzung des Trescaschen 
Kriteriums als Potential, doch erhalt 
man wesentlich andere Formen derh- und 
t-Kurven [vgl. 3 b)]. Die Geometrie des 
Vorganges |4Bt auf das plastische Po- 
tential schlieBen ! Auffalligerweise scheint 
h aber auch nach Mises eine Stelle star-- 
ker Kriimmung (die fast wie ein Knick 
wirkt) zu besitzen. 


Das Biegemoment M hingegen wird 
schon durch inkonsequente Benutzung 
des Trescaschen Kriteriums im Rahmen 
der Levy-Misesschen Gleichungen ver- 
ringert, hangt aber sehr stark vom Ver- 
festigungsgesetz ab. Die Abszisse der ~*y 025 05 O75 10 
Maxima ist variabel. Hierdurch unter- Aeneas 
scheidet sich wesentlich der Fall ebener — Abb. 5a—b. Spannungsverlaufe fir » = 2, 2 = 1 und Misessches Kriterium. 


tS 
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Formanderung (Proksa!) — wo die Maxima ziemlich unverandert bleiben — vom vorliegenden 
ebenen Spannungszustand. 

Natiirlich geht die Verfestigung unmittelbar in die Verteilung von Y ein. An ihr kann man 
(im Gegensatz zu M) sogar die Einfliisse von y und A getrennt beobachten. Leider ist Y ebensowenig 
wie o, oder dg direkt meSbar, so da® eine kiinftige Versuchsfolge ihr Hauptaugenmerk auf h, t 
und M zu richten hat. 


5. Elastische Riickfederung. a) Begriindung eines genaherten Ansatzes. Bei grofen 
plastischen Formanderungen gilt die Voraussetzung: 

Die durch das von uns errechnete Spannungssystem bedingte elastische Formanderung ist klein 
und andert die Querschnittsabmessungen kaum. 

Denkt man sich die elastische zur plastischen , 69 


Formanderung addiert, so geben sich nach (16) 3|% i 


nur geringe Anderungen der Spannungen und 
Krafte, die zu vernachlassigen sind. Auf diese 
Weise kénnen wir annehmen, die Spannungen 


seien richtig errechnet, und miissen zur erzwun- 2 
genen plastischen Formanderung y lediglich eine 


V8 


isa -4 
085 05 O15 40 0 025 YG ae 


Abb. 6. Spannungsverlauf fiir » = 2, 4 = 1 und Trescasches Kriterium. 


elastische Biegung y, hinzufiigen, welche durch das errechnete Moment bewirkt wird und nach 
Entfernung der 4uBeren Krafte zuriickfedert. 

Prinzipiell unzuganglich sind unseren Naherungen die durch Riickfederung im Werkstoff ent- 
stehenden Restspannungen. Sie bestehen wohl im wesentlichen aus Schubspannungen und ver- 
letzen unsere Voraussetzung 4. Zu ihrer korrekten Ermittlung miissen die Prandtl-Reufschen 
Seager sowohl fiir Biegung wie fiir Riickfederung zugrunde gelegt werden (vgl. Hill?, Ab- 
schn. ITI.2.). 


b) Formeln zur Berechnung der Riickfederung. Das Biegemoment M ist das durch 
unsere Rechnung vorgegebene Endmoment (26), y der (starrplastische) Biegezustand. Wir gehen 
zur Ermittlung von y, den iiblichen (auf der Bernoullischen Hypothese beruhenden) Weg (vgl. 
etwa Szabo* § 12; fiir eine genauere Analyse’, 10.6 oder Seames und Conway), miissen jedoch 


mene FuBnote 10 auf Seite 153. 

\ 2 Siehe FuBnote 3 auf Seite 153. 
3 Siehe FuGBnote 1 auf Seite 155. 
4 I. Szabé, Einfiihrung in die technische Mechanik, Berlin/Géttingen/Heidelbere 1958 
5 A. LE. Seames and H. D. Conway, J. Appl. Mech. 24 foe S. sa e 
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den gekriimmten Anfangszustand y > 0 beriicksichtigen. Dabei weicht die neutrale Faser von der 
Schwerpunktfaser ab. Statt dessen liefert die Bedingung der Kraftefreiheit bei linear verteilter 
Azimutalspannung 


; hx 


I—ya—*)* 
0 


A= (27) 


a 
f h 
) I-ya—5* 


Da y? < 1, ferner die ,,Wanderungsgeschwindigkeit“ dxy/dy <1, folgt bei Vernachlassigung von 
y? Oxy Oy 
M [l—y (l— xy)] 
Cua y : ie a am 
» f hx (x— xn) 
HU eaqa nk 


(28) 


Auf die etwas miihselige, aber triviale Ausfiihrung der Rechnung wurde verzichtet. 

Soll der bleibende Zustand y erreicht werden, so ist zunachst bis y + y, zu biegen, so dab sich y 
nach der Riickfederung einstellt. Die Gleichungen (27) und (28) sind i. a. numerisch auszuwerten. 
Sie gehen fiir y = 0 in bekannte Gleichungen iiber. 


6. Existenz und Eindeutigkeit der Losungen. a) Formulierung der Hauptergebnisse. 
Fiir |o,| << Y beweist man ohne Schwierigkeiten, da® die in 4 a) konstruierte Naherungslisung im 
Limes N-> co, Ay-> 0 ein Integral unserer Grundgleichungen darstellt. Sie ist die einzige, fiir 
welche A stetig und stetig differentiierbar von x abhangt, vorausgesetzt, Y hangt ebenso von A ab. 
Unter der Einschrankung |o,|< Y ist ferner Existenz und Eindeutigkeit der Loésungen des Tresca- 
schen Gleichungssystems gewahrleistet: Dieses fiihrte ja in Abschnitt 3 auf letzten Endes gewohn- 
liche Differentialgleichungen. Weitaus wichtiger scheint uns das folgende auf |o,| = Y beziigliche 
Ergebnis, welches zum AbschluB vorliegenden Aufsatzes bewiesen und diskutiert werden soll. 

Bei geeigneter Wahl von Y und h sowie hinreichend grofem y gestatten unsere Grundgleichungen 
keine Lésungen. 

Zum Beweis beschranken wir uns der Einfachheit halber auf das Trescasche Potential und geben 

= 1 vor, damit beide Seiten von (2’), (3’) zum gleichen Resultat fiihren. Ferner ecignet sich die stetige 


Funktion 


27 
Yen” (29) 
ee 8) 3S = r=.) 


(a = 1, b =5, p = konst.) als Beispiel — der Beweis kann ebensogut mit ahnlichen anderen, sogar 
mit stetig differentiierbaren Y gefiihrt werden, verliert dann jedoch an Einfachheit. 
Zum Vergleich mit 2 b) benutzen wir wieder unser urspriingliches Polarkoordinatensystem. 


b) Beweis bei Benutzung des Trescaschen Kriteriums. Zuniachst ist c, = ¢, [vgl. 2b)] 
unméglich. Denn o,, vom Rander = a her nach (3’) und (5) berechnet, liefert den Schnittpunkt 
i /3 mit der Kurve —Y(r) und erfiillt fiir nicht zu groBe r > /3 die Ungleichung o, < —Y. Da 
dies unmdglich ist, folgt ¢, < /3. Andererseits liefert o,, vom Rande r = b her errechnet, fiir hin- 
reichend groBe Konstanten p einen Schnittpunkt rs (3 < rg << 5) mit der Kurve —Y und erfiillt 
fiir nicht zu kleine r < rg ebenfalls die Ungleichung 0, << —Y. Es folgt ce, = rg > 3 > y3 = q. 

Aber auch c, > ¢, fihrt zum Widerspruch: 

1) Liegen die Stellen A = 0 in einem Teilintervall des Intervalles <c,, c.> dicht', so liegen auch 
die Stellen |o,| = Y dicht. Denn nach 2b) liegt zwischen je zwei Stellen r mit A = 0 eine Stelle 
\o,| = Y. Aus der Stetigkeit von o, und Y folgt o, = Y oder o, = —Y im ganzen Teilintervall — 
beispielsweise mége das erste gelten. An allen von r = 3 verschiedenen Stellen ist o, differentiier- 
bar: (2’), (3’) liefert den Widerspruch 

ar ¥) _ 


e os —2Y<—Y baw. Cp 8 prt SY. (30) 


1 Es mu dann unendlich viele geben — eine physikalisch untragbare Vorstellung. SchliefBt man sie von 
vornherein aus, wird der Beweis wesentlich einfacher. 
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2) Es enthalte nun jedes Teilintervall von <c,,¢,> ein anderes nicht ausgeartetes dy». das 
in welchem keine Stelle r mit A = 0 liegt. Das Intervall <d,,d,> sei maximal in dem Sinne, daB 


—O,75 


“65 


— 0,25 


0 025 05 OW 40 


0 0,25 05 0% 10 
Abb. 7b. 


Abb. 7a—b. Spannungsverliufe fiir » = 3, 4 = 2 und Misessches Kriterium. 


A=0 fir r=—d, und r=d,. Ferner 
kann <d,,d,> fremd zu r= 3 gewahlt 
werden: GiltA ¢ 0 auf dem ganzen Inter- 
vall <c,,c.> mit Ausnahme der Endpunkte, 
so darf dort nur eine einzige Form des 
Trescaschen Kriteriums richtig sein. 
Wegen (5) kame nur die dort nicht 
angefiihrte dritte Form |c,| = Y in Be- 
tracht, was zu (30) und einem Wider- 
spruch fiihrte. Es gibt somit eine Stelle 
r=caus dem Innern von <c,, ¢,> mit 
A=0. Dann wahlen wir <d,,d,> aus 
demjenigen der Intervalle <c,, ¢>, <c, ¢:>, 
welches r = 3 nicht enthalt. 

Auf <d,,d,> sind o, und Y iberall 
differentiierbar. Nach 2b) gibt es einen 
inneren Punkt de <d,, d,> mit |o,| = Y, 
z. B. o, = Y. Dann muBb ebenfalls 


@ eee 
Ge 

sein, weil sonst in unmittelbarer Nahe 
von r = d Punkte lagen, fiir die o, > Y 
wire. Jetzt aber kénnen wir (30) an- 
schreiben und einen endgiiltigen Wider- 


spruch feststellen. 


c) Verhaltnisse beim Mises- 
schen Kriterium. Ein entsprechen- 
der Beweis dirfte sich unter Benutzung 
des v. Misesschen Kriteriums formulieren 
lassen. Da er wahrscheinlich wesentlich 
verwickelter wird, unterdriicken wir ihn 
hier. Bemerkenswerterweise stellt sich 
der gewiinschte Widerspruch dabei 
nicht fiir |o,| = Y, sondern erst bei 


o,| =2 Y//3 ein, so daB ho, als Kurve 
ein legitimes Minimum [nach (16) etwa 
bei 6, = 0g = — Y| besitzen darf. 


d) AbschlieBende Diskussion. 
Der Beweis unserer Behauptung arbeitete 
unabhangig vom Verfestigungsgesetz. Es 
wird lediglich ein fester Biegezustand 
mit vorgegebenen h und Y betrachtet, 
wobei die Art der Entstehung gleichgiltig 
ist. Man kann y, h, Y z. B. durch span- 
gebende Formung und durch irgendwie 
inhomogen legierten Werkstoff erzwun- 
gen denken. 

In Wahrheit stellen sich bei konkre- 
ter Rechnung (so bei uns in Abschnitt 4, 
wenn y etwa denWert 0,9 erreichte. Wir 
geben daher nur Ergebnisse bis y = 0,7 
wieder) schon vor Auftreten des eigent- 
lichen Widerspruches groBe Ungenauig- 
keiten, negative Radikanden oder 
Doppeldeutigkeiten in z [wegen des 
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Minimums von o,; vgl. c)] ein, weil unser Beweis auch den allgemeinsten Fall mehrfacher Sprung- 
stellen von 0d ausschlieBen mufbte. Dieser ist aber praktischer Rechnung ohnehin kaum zuganglich. 


M 


ie tio % os 


O45 


040 


0,30|- 


0,25 


Abb. 8. Verlauf der Biegemomente. (Gestrichelte Kurve nach Tresca.) 


Der — bei Eindeutigkeit von z — wesentliche Widerspruch besteht darin, da nach dem FlieB- 
kriterium o, beschrankt bleibt, wahrend die Gleichgewichtsbedingungen (16) das Gegenteil er- 
zwingen wollen. Mit anderen Worten: Der wesentliche Widerspruch besteht zwischen FlieB- 


ti 


50 


40 


30 


60 


0 O25 05 075 10 
Abb. 9. Verteilungen von Y fiir y = 0,7. 
} 


| kriterium und Gleichgewichtsbedingungen; er ist unabhangig vom Stoffgesetz. Wir kommen zu 
- folgendem SchluB. 

| Die aufgetretene Unméglichkeit einer Integration unserer Grundgleichungen bedeutet, daB 
diese den plastischen Formanderungszustand in gewissen Bereichen falsch beschreiben. Griinde 
_ kénnten sein: 

| 1. Die Annahme ebenen Spannungszustandes, 

2. das Auftreten von Instabilitaten, 

| 3. die Vernachlassigung des elastischen Anteiles, 

4, ein prinzipieller Fehler der Theorie. 
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Die Griinde 1. bis 3. beziehen sich auf willkiirliche Vernachlassigungen im Ansatz. Fehlerquellen 
bei 1. und 3. sind unwahrscheinlich: 1. entspricht einer allgemein gebrauchlichen Annahme und 
kann durch ein Gedankenexperiment beliebig genau verwirklicht werden; 3. miiBte sich bei kleinen 
— nicht bei groBen — plastischen Forminderungen bemerkbar machen im Gegensatz zu den 
Erfahrungen bei unseren Rechnungen. Der Grund 2. kénnte verantwortlich sein. Allerdings 
miBte eine instabile Auslenkung des Balkens nach wie vor beziiglich O = 0 symmetrisch bleiben. 
In diesem Querschnitt miiBten also plétzlich starke Schubspannungen auf / wirken, die wir aber 
gedanklich durch Verkleinerung von h beliebig abschwichen kénnen. Unser Widerspruch dagegen 
ist unabhangig von h. 

Wir lenken die Aufmerksamkeit auf Grund 4. In der Tat wurde das FlieBkriterium bisher 
lediglich aus mathematischen Griinden gewahlt ohne Riicksicht auf andere bereits bestehende 
Spannungsbedingungen, so da® eine Unvertraglichkeit denkbar ist. Zur Klarung miissen grund- 
legende und quantitativ ecinwandfrei auszuwertende Versuche bei mehrdimensionaler Spannung 
und Formanderung dienen. 


(Eingegangen am 1, August 1958.) 
Anschrift des Verfassers: Dr. Horst Lippmann, Hannover-Doéhren, Am Lindenhofe 37, 
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Versagenszeit beim Kriechen 


Von F. Rimrott 


1. Einleitung. In dieser Abhandlung ist eine Beziehung (Kriechgleichung) zwischen dem 
TInnendruck in einer Hohlkugel und der Vergleichskriechgeschwindigkeit nach von Mises unter 
Verwendung von Ludwik-Dehnungen und -Spannungen entwickelt worden. 

In diese Kriechgleichung ist dann ein bekanntes, durch einen einfachen Zugversuch erhiltliches 
Kriechgesetz eingesetzt. Daraus ergibt sich eine Gleichung, die die Kriechgeschwindigkeit als 
Funktion der Kriechdehnung ausdriickt. Diese Gleichung kann integriert werden und gibt die 
Kriechdehnung in Abhangigkeit von der Zeit. 

Der letzte Abschnitt lehnt sich an eine von Hoff! entwickelte Theorie an, in der fiir den ein- 
fachen Zugstab gezeigt wird, daB die Zeitspanne bis zu einer unendlichen Kriechdehnung (Ver- 
sagenszeit) benutzt werden kann, um die Zeit bis zum tatsidchlichen Bruch (Bruchzeit) zu be- 
bestimmen. Die dickwandige Kugel kann in dhnlicher Weise wie der Zugstab behandelt werden. 
Nur miissen im Falle der Kugel (und des Hohlzylinders) zwei Umstande, namlich das Schwacher- 
werden der Wand und das Anwachsen der angreifenden Belastung, dadurch, daB sich die Innen- 
flache der Hohlkugel vergréBert bei gleichbleibendem Druck, beriicksichtigt werden, wahrend 
beim einfachen Zugstab nur Schwacherwerden des Stabes auftritt. 

Fur die in dieser Abhandlung vorgesehenen Kriechdehnungen und -zeiten kann Primarkriechen 
vernachlassigt werden. Das verwendete Kriechgesetz driickt das Verhalten von Werkstoffen 
wahrend Sekundarkriechens aus, und es ist auRerdem angenommen, da es das Verhalten von 
Werkstoffen wahrend Tertiarkriechens angemessen widergibt. 

Neben der ausfiihrlichen Behandlung der Hohlkugel, sind auch die Endergebnisse fiir Hohl- 
zylinder und Zugstab angefiihrt. 


2. Bezeichnungen. Wir benutzen folgende Bezeichnungen: 


Kriechbeizahl (Tag kg-” cm?"), 

Flache (cm?), 

Kriechmodul (kg cm~ Tag), 

Zuglast (kg), 

a’ Anfanglicher und verformter Innenhalbmesser (cm), 

AuBenhalbmesser (cm), 

h’ Wanddicke (cm), 

Kriechhochzahl, 

Innendruck (kg cm), ebenfalls anfangliche Zugspannung (kg cm~), 

r’ Halbmesser (cm), 

Zeit (Tag), 

Verlangerung (cm), 

Ludwik-Dehnung (Lingung), ebenfalls Vergleichsdehnung (Vergleichslangung), 
Kriechgeschwindigkeit, ebenfalls Vergleichskriechgeschwindigkeit (Tag), 
Ludwik-Spannung (wahre Spannung), ebenfalls Vergleichsspannung (kg cm~*). 


Q Or eye SS SSS eye aa 
(ep 


Die FuBzeiger bedeuten: 


am Innenhalbmesser, 

am AuBenhalbmesser, 

radial, 

beim Versagen, 

tangential (Kugelkoordinate 7), 
tangential (Kugelkoordinate @). 


1 N. J. Hoff, Appl. Mech., Trans. ASME 75 (1953), S. 105. 
12 
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3. Annahmen. Die spater folgenden Ableitungen sind auf folgende Annahmen aufgebaut: 
1) Die Wichte des Werkstoffes bleibt unveradndert: 
Ey +é& +&, = 0. (1) | 
2) Die Verhaltnisse der Hauptschubspannungen zu den Hauptschubdehnungsgeschwindigkeiten 
sind gleich: 


06—Or _0r—Oy _ Sy— 98 (2) 
éo— & &&—eEy éey— és’ 

3) Es gibt eine Beziehung zwischen der Vergleichsdehnungsgeschwindigkeit und der Vergleichs- 
spannung, die gleich der Beziehung zwischen Dehnungsgeschwindigkeit und Spannung beim ein- 
fachen Zugversuch ist: 

é=f(o). (3) | 


4) Der Kérper ist aus einem isotropen und homogenen Werkstoff hergestellt und andert seine _ 
Gestalt wahrend des Kriechens nicht. 

Ahnliche Annahmen wurden auch von Soderbergi, Marin?, Bailey®, Johnson* und Nadai? | 
gemacht, mit dem Unterschied, da8 in dieser Abhandlung Ludwik-Spannungen und -Dehnungen 
benutzt werden. 

Die Annahmen 1) und 2) kénnen zusammengefaBt werden zu 


: 1 1 : 1 
Eg = 7 |% — z= , + oy) ’ eee 


1 


c—Hy tal), yoke bere), a 


wobei der Kriechmodul K durch 
K=< (5) 


é 
gegeben ist und unter Benutzung von Annahme 3.) fiir jeden Spannungszustand berechnet werden 
kann. 


4, Kriechgleichung. a) Dickwandige Kugel. Der Kiirze halber sei die Beziehung zwischen 
dem Innendruck und der Kriechgeschwindigkeit hier Kriechgleichung genannt. Die Spannungen, 
Dehnungen, Kriechgeschwindigkeiten sind einerseits voneinander abhangig, namlich o = f(é), 


und ¢ = f(é), andererseits sind alle drei GréBen auch Funktionen vom Halbmesser r und der Zeit t, 
dae == f(r, t). ast. 


Lp Antangszustand — Verformter Zustand 


Abb. 1. Element einer verformten Kugel. Abb. 2. Querschnitt einer dickwandigen Hohlkugel 
oder eines dickwandigen Hohlzylinders. 


Die Bedingung fiir radiales Gleichgewicht einer Kugel unter Innendruck (Abb. 1) ist fir den 
verformten Zustand gegeben durch | 


, do, 


rs, =2(%—<4,). (6) 
1 C. R. Soderberg, Trans. ASME 58 (1936), S. A 150. 
J. Marin, Trans. ASME 58 (1936), S. A 21. 
3 R. W. Bailey, Proc. Institution Mechanical Engineers, England, 154 (1951), S. 425. 
. s E. Johnson, Proc. Institution Mechanical Engineers, England, 164 (1951), S. 432. 


- Nadai, J. Appl. Physics 8 (1937), S. 418. 
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Die Ludwik-Dehnungen sind (Abb. 2) 


é =In (7) = In(1 ue = (7) 

Rs In] an (1 a Z| (8) 

= m(*) =In (1 a a (9) 
Die GréBe u 148t sich aus den Gleichungen (7) und (8) entfernen, wenn diese zusammengefabt 
werden zu 

r $2 — exp (c, — 65) —1. (10) 


Die Beziehungen in Gleichungen (7) und (8) kénnen benutzt werden um die Gleichgewichts- 
gleichung (6) umzuwandeln zu 

00; 

"Or 


= 2 (0) —o,) exp (e, — 6). (11) 


Eine hohle Kugel unter Innendruck gehért zu den Fallen, in denen die weitere Ableitung der 
Kriechgleichung entweder auf die Vergleichskriechdehnung oder die Vergleichskriechgeschwindig- 
keit aufgebaut werden kann. Die Vergleichskriechdehnung nach von Mises, die hier benutzt werden 


soll, ist festgelegt als 


ce ie (Eo é,)* + (e, Sees) Shey es) (12) 


Mit Hilfe der Gleichungen (7), (8) und (9) sowie unter Beriicksichtigung der Bedingung (1) laBt 
sich (12) vereinfachen zu 


&E=2& =—2e,=—é,. (13) 
Die Vergleichsspannung nach von Mises ist 
1 = 
i iz /(os —a,)? + (, Oa Og Oo) a (14) 


Da nach (7) und (9) ¢, = é, ist, ergibt sich aus (4), daB og =o, ist; damit kann (14) vereinfacht 
werden zu 
OG 9 ae On Oy) ae Orn (15) 


Die Einfiihrung der GréBen o, gegeben durch (15), und ¢, gegeben durch (13), in (11) fiihrt zu 
oe 2 exp ic | ; (16) 


r aT 
or 


Gleichung (10) kann nach Einfiihrung der GréBe ¢ aus (13) integriert werden und ergibt 
2 T(t) 
e=—ln ( fe ra (17) 


wo T(t), die Integrationskonstante nach der Veranderlichen r, eine Funktion der Zeit allein ist. 
Fiir die Vergleichsdehnung «, fiir r = a, wird (17) 


e, = In (1 a | : (18) 


a? 


Durch Zusammenfassung der Gleichungen (17) und (18) kann man T(¢) ausscheiden und erhalt 
a\3 Bea 


eine Gleichung, die die Verteilung der Vergleichsdehnung tiber die Kugelwand fir eine gegebene 


Vergleichsdehnung ¢, am inneren Halbmesser darstellt. 
Die Kriechgeschwindigkeit é ist die partielle Ableitung der Dehnung ¢ nach der Zeit: 


ane (3 €a 
(F] eis (| ee (20) 


Z 
S Sys 


Na 
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Die partielle Ableitung 9é/gr ergibt sich aus (20) bei Beriicksichtigung der Beziehung (19) zu 


0€ one 3 ¢€ 
bg se eyes 21 
soa esl ep 
Aus (16) und (21) erhalt man 
do, 2 6 
dey ane co 


nach Integration zwischen entsprechenden Grenzen ergibt sich die Kriechgleichung fiir die dick- 
wandige Kugel zu 


Ea 


p= | +d, (23) 
& 
wobei 
ane 3 Eq) 
Bi exp ( 5) } 
wes (24) 


ist, wie sich aus (20) ergibt. 


b) Dickwandiger Zylinder. Fir einen dickwandigen Zylinder mit geschlossenen Enden 
unter Innendruck kann die der Gleichung (23) entsprechende Kriechgleichung auf ahnliche Weise? 
abgeleitet werden. Sie ist 


Daa < di (25) 
mit 


= 


1+ (2) fe (5 «—1) oa 


c) Dimnwandige Kugel. Die Kriechgleichung fiir eine diinnwandige Kugel laBt sich leicht 
ableiten. Wenn die Wichte des Werkstoffes unverandert bleibt, dann muB sein 


Aga*h’ =4na%h. (27) 
Aus (7) und (13) ist ersichtlich, daB 


a =aexp— (28) | 


ist. Demnach wird gemaB (27) und (28) 
h’ = hexp (—e). (29) : 


Fir eine diinnwandige Kugel unter Innendruck ist angenahert | 


mi | 
ap } 
Og = gape ae Oy (30) 
Anfangszustand — Verformter Zustand 
und 
Abb. 3. Querschnitt einer diinnwandigen | 
Hohlkugel oder eines diinnwandigen Zylinders. Of, = Oe (31) | 


Aus (15), (28), (29), (30) und (31) ergibt sich die Kriechgleichung fiir die diinnwandige Kugel zu | 
2h 3¢€ | 
| (32) | 


d) Diinnwandiger Zylinder, Die Kriechgleichung fiir einen diinnwandigen Zylinder mit | 
geschlossenen Enden unter Innendruck 1a8t sich in ahnlicher Weise wie Gleichung (32) ableiten. | 


1 F, Rimrott, ASME-paper 59—APM-6. 
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_ Fiir den diinnwandigen Zylinder ist 


pS a a exp (_ 3s), (33) 


e) Zugstab. Fir einen Stab unter konstanter Zuglast! laBt sich die Kriechgleichung ableiten zu 
p=o e-= (34) 
mit p = P/F, wo P die konstante Last und F, der Anfangsquerschnitt ist. 


5. Kriechgesetz. Das am meisten benutzte Gesetz fiir die Beschreibung des Kriechverhaltens 
_ eines Werkstoffes lautet 


é = Ao’, (35) 


/ was nun in die in den vorhergehenden Abschnitten abgeleiteten Kriechgleichungen eingesetzt 
_ werden wird. 


a) Dickwandige Kugel. Wenn (35) in (23) eingesetzt wird, dann ergibt sich nach Aus- 
 fihrung der Integration 


Zier : 
ee 3 All2 (oP =S 3/") . (36) 


Driickt man ¢é, durch (24) aus, so kann (36) nach é, aufgelést werden: 


(sao) YT" 
b 25 
eS ee : (37) 
Ms 2”? nt a \8 3 &@ 
(GY fee) 


Gleichung (37) ist eine Differentialgleichung, die gelést werden muB, wenn man eine Beziehung 
zwischen Kriechdehnung und Zeit erhalten will. Integration ergibt fiir konstanten Innendruck 


€ 
a al an ( 2 1/n7]n 
i e 3 fipene Ralry (38) 


alee | 


a\3 oe ( = - 
b Boe 
. de (39) 


Mit der Abkiirzung 


wird die Zeit 
n n 
pees (40) 


Das Integral J (39) kann auf folgende Weise ermittelt werden. Setzt man 


ons “i 


dann wird 
Ye 
fea 8 le Li (42) 
3 (Le yyy 
af 
wobei 
1 = (a/b)8” (43) 


1 Siehe FuBnote 1 von Seite 169. 


| 
| 
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i: 3 Ps E 7 1/n | 
— e oe } 
a aaa EEE (44) | 


der Grenze 0 und 


der Grenze ¢, in Gleichung (39) entspricht. pe 
Beschrankt man die Werte fiir n auf ganze positive Zahlen, dann lat sich das Integral in die 
folgende binomische Reihe entwickeln: 


n—1 


Ya 
Baa eres laine) eos 45) |} 
Abeer Z eee a Co) | 


Pik § 


Das Integral kann nun Glied fiir Glied berechnet werden. Jedes Glied ist von der Art 


Jo=—(4) | ae Wo 


und ergibt? 


Jn = —( : (rine 1) fee TG + 1) 


m-+ 1 
(ea lYe 
= =. (R cos eee — S, sin eee || (47) 
vy=1 ly, 
mit n =20—S, Si— 0 nod cries Sov ] 
ley hstiti=s,0, 1,52. WS eds gotta k 
und 
il 22 
Razin yt —2y coe 4.1), (48) 
2vn 
y— cos S 
TORE Ue IE 
2van ; 
sin {—— 
Ce) | 
Fir m = — 1 wird das Integral (46) 
fg eee (49) 
tae a be mi 


Wenn man die Integrationen nach Gleichungen (47) und (48) so oft wie nétig durchfithrt und 
die Ergebnisse addiert, so ergibt sich 


al 


_ Sip le yee ul % n 
Fah vit ees 


m=0 
2S 2v(m + 1)z . 2v(m+1)a 
+ ES ante aee ml 


wobein > 2 und eine ganze Zahl, y, y, und y. durch (41), (43) und (44), m, 0, s, R, und S, durch (48) 
gegeben sind. Fiir den Fall, daB n < 2 ist, 1aBt sich das Integral (42) auf andere, in jeder Integral- 
tafel auffindbare Formeln auswerten. Der aus (50) erhaltene Wert fiir J muB in (40) eingesetzt — 
werden, um eine Beziehung zwischen Kriechdehnung und Zeit zu erhalten. 


=In(y—1) + 9 ay +) 


= (50) 


v1 | 


b) Dickwandiger Zylinder. Fiir den dickwandigen Zylinder erhalt man an Stelle von 
Gleichung (40) 


n 
pee, 
y3" 4 p® 
1 W. Grébner u. N. Hofreiter, Integraltafel, Erster Teil, Seite 18. Wien und Innsbruck 1949, 


(51) 
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[ 
he 
h 
| 


und 


i! 
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maak re (1 — x)" dx 
a\2/n 
St Pm (5) (53) 
el exp (3 &0) 
(54) 


eV =i] | 


a) 


Das Integral (51) 1aBt sich ahnlich wie Integral (42) lésen. 


|- c) Dinnwandige Kugel. Fiir die dinnwandige Kugel ergibt sich die Zeit als Funktion der 


Vergleichsdehnung zu 


(55) 


d) Dinnwandiger Zylinder. Fiir den diinnwandigen Zylinder ergibt sich folgende Gleichung: 


exp (= \3n e) 


2h 


y3A el 


e) Zugstab. Fir den Zugstab ist! 


as 1 — exp (— ne) 
a n A p® ; 


(56) 


(57) 


Die Gleichungen (55), (56) und (57) ergeben sich als Lésung der Differentialgleichungen (32), 


(33) und (34), die entstehen, wenn o durch A—!/" ¢1/" aus (35) ersetzt wird. 


6. Versagenszeit. Die Versagenszeit t, sei festgelegt als die Zeitspanne zwischen dem Beginn 
des Kriechens (bei dem ¢ = 0 und t = 0) und dem Erreichen einer unendlichen Formanderung 
(gegeben durch ¢ = co und t = ,). 


Vergleichsdehnung &, amInnenhalbmesser 


h/a=1 oder bla=2 


"GF, 


Abb. 4. Verformung einer dickwandigen Kugel in Abhingigkeit von der Zeit. 


° 


50/0 
der Versagenszeit ty nach Gleichung (58) 


a) Dickwandige Kugel. Die Beziehung zwischen der Vergleichsdehnung ¢, am Innenhalb- 
messer und der Zeit ist in (40) gegeben. In Abb. 4 sind Kurven — Vergleichsdehnung ¢, gegen 
Zeit t — aufgetragen. Aus ihnen ist ersichtlich, daf die Zeit, in welcher Bruch zu erwarten ist 


1 Siehe FuBnote 1 von Seite 169. 
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(Bruchzeit) und die Zeit, fiir die ¢, = © ist (Versagenszeit t,) so nahe beieinanderliegen, daB die | 
Zeit zum Bruch geniigend genau durch die Versagenszeit angenahert werden seat Nimmt man — 
zum Beispiel an, da® Bruch eintritt wenn ¢, = 1,0 ist, dann ist die Bruchzeit 54% der Versagenszelt ' 
fir n = 1, 77% fir n = 2, und 94% fir n = 4. Je groéBer n, desto besser ist die Annaherung. 


Tage) 1=4 
A=1- 10" cmb kg * Tag? 


Theorie fir 
diinnwandige Kugel 
dickwandige Kugel 


10° 


te 


Versagenszeit ty 
7 
Ss 


10 


ONG oF Q1 10 70 


Wanddicke/Innenhalbmesser ha 


Abb. 5. Versagenszeit fiir Kugel unter Innendruck. 


Wenn man nun bedenkt, daB die iiblichen Werte fiir n im Kriechgesetz (35) 4 und gréBer sind, dann 
ist der Fehler bei der Benutzung der Versagenszeit zur Bestimmung der Bruchzeit, fiir den F all, 
daB die Ludwik-Dehnung beim Bruch ¢, = 1.0 ist, ungefahr 6% oder weniger. Die Versagenszeit 


kann also gut fiir die Bestimmung der Bruchzeit verwendet werden. Fiir sie ergibt sich aus (40) 
fiir die Versagenszeit 


20 nt 


= 5a gone (58) 
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wobei 
co 


1( ee 3 &4 1/n7]n 
dala ale 


i= ere ae (59) 
THe leo el 
6 
was nach Einfiihren der Veranderlichen y (41), die folgende Form annimmt: 
1 
2n [ (—y)" dy 
ah — ae = 
3 | Ce ean ae ea 2 eo) 


(a/b)3/" 


Kurven fiir die Versagenszeit t, in Abhangigkeit vom Wanddickenverhaltnis h/a sind in Abb. 5 
gezeigt fiir die dort angegebenen Kriechkonstanten. 


b) Dickwandiger Zylinder. Die Versagenszeit fiir den dickwandigen Zylinder ist gemaB (51) 


n® 
b= aap” es 
mit 
1 
Bee (1 — x)” dx 
ee Ue be 


(a/b)?! 
c) Dinnwandige Kugel. Fiir die diinnwandige Kugel ergibt sich die Versagenszeit aus (55) 
indem man € = © setzt, zu 


22h 
°"3nA (ap) ° 


(63) 


Gleichung (63) ist ebenfalls in Abb. 5 graphisch dargestellt. 


d) Diinnwandiger Zylinder. Fiir den diinnwandigen Zylinder erhalt man die Versagenszeit 
indem man ¢ = 0 in Gleichung (56) einsetzt: 
a ebeeh ane oe 
asm A (V3 ap) 2 


e) Zugstab. Die Gleichung fiir die Versagenszeit eines Zugstabes unter konstanter Last wird 


zufolge (57) 


(64) 


1 
ang pie (65) 

7. SchluBbemerkung. Eine Hohlkugel unter Innendruck wird, wie bereits erwahnt, natiirlich 
etwas eher bersten als die Versagenszeit angibt. Bersten tritt ein, sobald die Wanddicke so klein 
und die Belastung so groB ist, daB die Spannung in dem Werkstoff diesen zum Bruch bringt. 

Die Versagenszeit ist auch dadurch beeinfluBt, daB sich die Kugel infolge unvermeidlicher 
Umstande, wie Ausmittigkeit, ungleiche Wanddicke, Werkstoffehler, Inhomogenitat und Aniso- 
tropie, nicht symmetrisch ausdehnen wird, wie das in der Theorie angenommen ist. Die soeben 
erwahnten Fehlerquellen treten jedoch erst dann merklich auf, wenn die Verformung dem Ende 
zugeht und zu einer Zeit die der Versagenszzit so nahe liegt, daB ihr KinfluB auf das Endergebnis 
als unbedeutend angenommen werden kann. 


8. Zusammenfassung. Durch die Verwendung von Ludwik-Dehnungen und -Spannungen ist 
es in einigen Fallen méglich, eine Bestimmung der Zeit bis zum Bruch eines kriechenden Werk- 
stiickes auszufiihren. Zu den Fallen, in denen sich eine solche Bestimmung durchfihren laBbt, 
gehéren Hohlkugel und Hohlzylinder unter Innendruck sowie der einfache Zugstab. In dem vor- 
liegenden Aufsatz ist der Fall einer hohlen, dickwandigen Kugel unter Innendruck ausfiihrlich be- 
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handelt und sind die Endergebnisse fiir den dickwandigen Zylinder, die diinnwandige Kugel und | 


den diinnwandigen Zylinder sowie den Zugstab kurz angefihrt. 
Der Grundgedanke der hier entwickelten Theorie ist, die Dehnung in den Kriechgleichungen 


nach unendlich gehen zu lassen. Die Zeitspanne, in der eine unendliche Verformung erreicht wird, | 


ist Versagenszeit genannt worden. Es wird gezeigt, daB die Verformung am Ende so ungeheuer 


schnell vor sich geht, daB die Zeitspanne, um enendliche Dehnung zu erreichen (Versagenszeit), und | 
die Zeit bis zum tatsachlichen Bruch (Bruchzeit) nur geringfiigig voneinander abweichen, und daB — 


daher die Versagenszeit zur Bestimmung der Bruchzeit verwendet werden kann. 


(Eingegangen am 14, August 1958.) 


Anschrift des Verfassers: Prof. Dr. Ing. Friedrich Rimrott, University Park, Pa. (USA), Pennsylvania State | 


University, Abt. fiir Technische Mechanik (Prof. Dr. Marin). 
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Thermal stresses in a heat-generating cylinder: 
The variational solution of a boundary layer problem 
in three-dimensional elasticity 


By G. Horvay, I. Giaever, and J. A. Mirabal 


1. Introduction. When self-equilibrating radially symmetric normal and shear tractions s,(r), 
t,(r) act on an end-face of a cylinder, the ensuing stresses and displacements are concentrated in a 
thin boundary layer of the cylinder, about two radii deep. The stress distribution may be calcu- 
lated by using a pair of Sadowsky-Sternberg axisymmetric stress functions _,(z, r), D,(z, 7) 
which identically satisfy the equilibrium equations. We write ,, D, as products, g,(z) f,(r), 
G,(2) F,{r), where f,(r), F,(r) are polynomials of degree k + 2 (k = 2, 4, 6,...) which satisfy certain 
boundary conditions and orthogonality relations ; their appropriate derivatives constitute the two com- 
plete sets of end tractions s,(r) and t,(r). The axial functions g,(z), G,(z) — decaying with distance from 
the loaded end as e~%* *, where a, is the real part of an eigenvalue parameter y, — are determined 
from Euler-Lagrange equations. It is found that the variational approach used in this fashion leads 
to displacements in the fundamental mode which satisfy the compatibility requirements within, 
almost, graph paper accuracy. The theory developed is applied to determining the stresses in a 
cylinder of unit radius and axial length L greater than five radii, when a radial temperature distri- 
bution 0 = @,, (1 — r?) is maintained. The surprising result is obtained that the maximum stress 
in the cylinder, a longitudinal stress o, = + 0.501 Ea Oy/(1 — ) (plus at r = 1, minus at r = 0), 
develops at a relatively large distance, 2.3 radius, away from an end-face. This barely exceeds 
the value + 0.500 E « Oy/(1—¥v) which prevails throughout the remainder of the cylinder 
2.9525 L—2.5. 


2. The Rigorous Solution, using potential Functions. We are considering the problem of deter- 
mining the stress distribution in an elastic cylinder which is free from tractions on its curved surface: 


Tz. 1) == () 5 (la) 
o,(%, 1) = 0 (1b) 


and is subject to radially symmetric shear end tractions t°(r) and normal end tractions o;(r) 
which are self-equilibrating, i. e., which satisfy 


fo rdr=0. (1c) 


We refer to this problem as the “end problem’. The symbols o,, o,, gg denote the axial, the radial 
and the hoop stress components, 7 is the zr shear stress, and w and u the axial and radial displace- 
ments (see Fig. 1). We use the cylinder =. 

radius as the distance unit, the double 
shear modulus as the stress unit. We use 
the word ‘“‘traction”’ in the sense of “boun- 
dary stress’. The superscript ° refers to 


0 
the z = 0 cross section. If, as we assume, aan 
the length of the cylinder exceeds four, 
the interaction of the effects at the two - 
ends may be neglected, and the cylinder Nes 


may be regarded assemi-infinite. 
In a previous paper? — to be referred Pog Geyer 
to hereafter as 1 — it was shown that the 


Sadowsky-Sternberg stress functions ¢(z, r), é 
@(z,r) are ideally suited for solution of the axisymmetric “end problem.’ Indeed, on 


1 G. Horvay and J. A. Mirabal, The End Problem of Cylinders, J. Appl. Mechanics 25 (1958), p. 561. 

2A he poe eaenal approach is to use a “‘Love” function (see H. Poritsky and M. F. Bolster, Thermal 
Stresses in Cylinders, GE Report GE-DF-52 GL 128, and Refs. 3—5 of I). While, naturally, both methods 
lead to the same rigorous solution problem, only the Sadowsky—Sternberg functions are adaptable to a variati- 


onal approach. 
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writing 


Ope 2 Os ata Lee) O,, +17 ®,,; 


1 
65— V7 @ me f (1 2) ®.,, 


(2) 
0) SA at ed et ee rete eel | 


tT=—,,+7®,,,—2(1—»)@,,—rV’®, J 
it is found that the equilibrium equations 
= x ! Gr — 8. ae 0 : | 
ia 3B | 
eee er (3) | 
(e)-4 F or + r 


are identically satisfied!. One also notes that the stress Expressions (2) do not depend expli- 
citly? on z. On writing the product function representations 


y(z, 7) = g(z) f(r)» D(z, r) = G(z) Fir) . (4) 
and introducing the notation 
> lade 
/ — se 4p = # ete., (5a) 
trr=f', tp=rF’4+(3—2»)F’, (5b) 
/ t , tr 
ss= | f Sp= tpt (5c) 
the stress expressions become 
o,=fe" +if'g +r F + (1—29) FIC”, (62) 
Op =fe- +f’g+(1—2») Fe", (6b) 
O, = srg t+spG, (6c) 
t=—tg—t,C’. (6d) 


The next objective is to determine g and @ so that the stress-strain relations (the “‘compati- 
bility”’ relations) 


(+) =6,—v(+0,), (7a). 
(1 +») — =o9—» (6, +4,), (7b) 
(1 +) =0,—9(6, +07) (70) 
4 Sars ea 


be also satisfied. This is accomplished, as is well known, by choosing the functions gw, ® as har- 
monic functions 


VYp=0, VWe=0 (8) 


* In their original form, as given by Sadowsky-Sternberg (Lecture Notes on Three Dimensional Elasticity, 
by M.A. Sadowski, unpublished), the 7? terms were absent from (2), and hence (3) was satisfied only on 
condition that (8) hold. The \/* terms necessary to insure satisfaction of (3) without the requirements (8) 
were first introduced in 1. 

® This as will be shown, insures an exponential decay of the stresses with z. The more familiar Boussinesq 
representation of the stresses involves z explicitly. It is therefore unsuitable for the treatment of the “end 
problem”, even after appropriate \/* terms have been incorporated in the stress expressions, so that (3) be 
satisfied. However, the Boussinesq representation is the proper one to use if the lateral surface of the cylinder 
is loaded by self-equilibrating tractions o,(z, 1), t(z, 1) and the end-faces z = 0, L are free. 
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and defining the displacements by 
UO PO, , w= —p@—r1O,, = 4(1 =) @,. (9) 


We are interested in exponentially decaying solutions! of (3) and (7), i. e., of (8). Such are g, 
@® where 


ge)=E@) =", f=, FO) =BIg(yr). (10) 
- It is seen that the ensuing stresses o,, 09, o, are even function in r, tT is an odd function. The choice 
ie ee EO RAD oe 
insures that the boundary conditions (1a, c) are satisfied?, and the choice of y as the root of 
1c irs (2) 
insures that the condition (1b) are also satisfied. We denote the roots 
Ve =O, + 1B, (13a) 
which have positive real and imaginary parts by y2, y4,... Poritsky-found that? 
Yo = 2.72212 + 1 1.36209 , (13b) 
VY, = 6.05962 + 7 1.63791 . (13 e) 


It is seen that —y, and + y,, conj are also roots. The functions y and @ associated with y, are 
denoted by @,, ®,. Because y, is a complex number, g, and ®, are complex functions. Their 
real and imaginary parts (i and §) individually satisfy (3) and (7). 


3. The variational method. Through Expressions (10)—(13) we have acquired two infinite 
families of stress function pairs which identically satisfy Eqs. (3) and (7). These are the real parts 
of gy, and ®,, and the imaginary parts of y, and ®,. These functions may be adapted to solution 
of the end problem: the determination of the stress distribution in the cylinder under given self- 
equilibrating end tractions. One only has to take suitable linear combinations of the stress function 
pairs, so that their appropriate derivatives reduce to (or approximate) the given end tractions, 
and thenone obtains an exact (or an approximate) solution to the end problems. 

But the use of the functions Jig,, iO, Sy, GO, encounters great practical difficulties. 
The Besse functions of complex arguments are inadequately tabulated. Furthermore, the func- 
tions y,, Dl, are not bonafide product functions: their product nature disappears on taking real 
and imaginary parts. Because of this, the functions give rise to both o;(r) and t°(r) boundary 
stresses, which are hopelessly intermingled: ordinary Fourier expansion methods for representing 
given boundary tractions in terms of characteristic tractions s,(r), t,(r) do not apply. 

By sacrificing the requirement (8) for harmonicity we may salvage both the bonafide product 
representation (4) of the stress functions, and the method of individual Fourier expansion for each 
end traction o,(r) and t°(r). This is accomplished by the use of the variational method, to which 
we now proceed. The functions y,, ®, we shall define presently automatically satisfy the equili- 
brium equations (3), but only approximately the compatibility equations (7). Indeed, the quality 
of the variational approximation may be judged by the degree of adherence to the 
compatibility conditions (7). When the variational stress functions @,, ®, are introduced 
into (7), and the radial displacements u, as determined from (7a) and (7b), respectively, are found 
to be barely distinguishable, and the axial displacements w determined from (7c) and (7d), respec- 
tively, are likewise found to be barely distinguishable, then we may regard our solution ~,, ®, an 
excellent approximation. If the agreement is only fair in the displacements, the average of the two 
displacement determinations may still be quite representative, and the stresses may still be quite 
reliable. If the agreement of the displacements, calculated by the two procedures, is poor, then the 
stress distribution may be judged as having only qualitative usefulness. 


1 There exist also solutions with e+ 7# variation. They enter the analysis only in cases where the cylinder 
is so short (L < 4) that the interaction of the effects from the two end-faces cannot be ignored. 


2 Tn fact, etc. f or) r dr = 0 is satisfied for all z. 
0 


3 All numerical results stated in this report are based on a value of the Poisson ratio y = 0.3. 


0 (Or) 4 0;(0)&t(Or) 
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We shall show that the approximation of 1, based on the assumption 


Method A: either o== 0; D = G(z) F(r) (14a) 
or yp =s8la)f(r), P=0 
belongs into the last category, the approximation based on the assumption 
Method B: yw = G(z) f(r), ® = G(z) Fir) (14b) 
(i. e., g =G) belongs (for the fundamental traction) into the intermediate category, and 
Method C: gp =g(z2) f(r), D = G(z) F(r) (14:c) 


(with the g, G relation established by the calculus of variations) belongs into the first category for 
the fundamental normal traction, and into the second category for the fundamental shear traction. 

As an illustration — and in anticipation of the results to be arrived at in this report — we 
compare in Fig. 2(a) the variationally determined end-face hoop stress o@(r), produced by the 
“fundamental normal traction mode” 


Off) = 2rt—1  r (r) 0 (15) 
with the rigorous solution of (2), (7) which Poritsky obtained when he used a slightly modified set of 


boundry aconditions.t The boundary tractions (15) and Poritsky’s modified tractions are illustrated 
in Fig. 2(b). It is clearly difficult to judge whether the slight discrepancy between the o4(r) 


0 


| if 
T 05 10 


A: 9 =0,9-F(r)G(z) /, HA 0 
B: p =F (r)G(Z), 9=F(r)G(z) ~ { 
C: p =f (r)g(z), =F (r)G(z) 
P: Poritsky ~ al 


~ 


,(07)86 (OF) 


~__ ae 


Vig. 2. The hoop stress o9(0,r) calculated by the variational Methods A, B, C for the rigorous boundary conditions 
oz (0,r) = 2r*— 1, x(0,r) = 0, is compared, Diagram (a), with rigorous (Poritsky) solution for the approximate boundary 
conditions illustrated by dashed lines in Diagram (b). 


curves C and P should be attributed to the approximate nature of solution C or to the approximate 
nature of boundary conditions P. : 

Our variational approach involves the following steps. We introduce the functions (4) into the 
complementary energy expression 


oo 1 
Us 2m fds [lag py + 06 + o,)?—9, 069 — (0, + 69) 6, + t2|rdr (16) 
0 0 


1 In his calculations Poritsky employed a biharmonic stress function of form 

G@=(AR+ BRL, +(CR4+D9)L, 
where L,, Ly are the Love functions belonging to Yoo Ya, and A, B, C, D are coefficients determined by least 
square fitting of the boundary stresses to the tractions (15). 
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and, regarding f(r) and F(r) [or, more properly, f,(r) and F;,(r)] as prescribed quantities, we integrate 
out over the variable r. That gives us 


1+ yy 7 , vt 0 ’1 
ga UT = JX (6 8',8"s 66,6") dz, (17a) 


& = doo 8? + a1g87%+ Qg2g 7 + Qo. gg’ + Co G? + ¢4, G2 + Coo G''* + Oyo GG” 
4 + boo gG + big GC’ + bys 8” Ge + bag G" + boy g’ G (17b) 
where 
1 


Ao = (l—») f sfrdr, 


0 
i 


a,=—(1+7)f#rdr, 


0 


: (18a) 
a, =(1—») ff?rdr, 
0 
1 
Ay, = —(1— 3») f#rdr, 
6 
1 
0 
1 
6 
1 
bo, = (1—») [flr F’ +2(1—2») F]rdr, (18b) 
1 1 
bog = 24 ftetprdr+2y(l—») ft; F’rdr, 
3 0 
1 
bog =20f tptprdr. 
0 
1 
¢ =z | 4 d 
00 9 nT Ok 
0 1 
i= (L+9) [thrdr, 
ey i <0 1 (18c) 
cay [FPrtdr + (1-29) 1 —») | FUR r+ 12) Fra, 
0 ) 
1 1 
con = 9 | thrdr—208 | te Fr dr. 
3 0 


One notes that the functional & in (17a) does not involve z explicitly, so the Euler equations derived 
from (17), whether we place g = 0 or G= 0 as in Method A, or write g = G as in Method B, or 
treat g(z) and G(z) independently, as in Method C will be differential equations with constant 
coefficients. This signifies an exponential variation for the axial functions. 

Of great importance in the variational procedure is the judicial choice of the radial polynom- 
ials f,(r), F,(r) for which the g,(z), G,(z) functions are to be determined; or, what amounts to the 
same thing, the proper choice of the characteristic boundary tractions s,(r), t,(r) for which the 
stress functions ¢,(z, r), D,(z, r) are to be determined; or, what amounts to the same thing, the 
proper definition of the rule of orthogonalization for the polynomials f,(r), F,(7). 


4. Method A. The approximation (14a) was carried out in 1. In that paper we recognized that 
the lowest degree normal traction polynomical (non identically zero) which obeys (1c) is 


s,(r) =272—1. (19 a) 
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The subscript k = 2 refers to the degree of the s(r) polynomial. We use bold face letters s, Peele 
to indicate that the normalization 
GW sil (20) 


is used. Integration of (5c) then gives 


f(r) = (19b) 


and (19b) is seen, by (6d), to obey (la). (The integration constant is zero since the shear stress 
is an odd function of r.) A further integration of (5b) now leads to 


f,(r) = (1— 8)? (21) 
for the case ® = 0, and to 
8—9rv+27? We r4 
Fett) = 76 (1) Oe Oe a= Be, 


or 


F,(r) = g5.9 [r! — 3.176471 1? + 8.058823] (22b) 


for the case p = 0. [The integration constant is determined by the boundary condition (1b).] 
Introducing (19), (21), (22) into the strain energy expression (17) there results an Euler diffe- 
rential equation of the form! 


WeY— Be’ +Cg=0. (23) 
Its solution may be written as 
for o,(r) = given , c(h) 0; pe) ese" (cos pat ai sin | ; (24a) 
for rr) =given, o'r) = 0: ——(z) =e ware (24b) 
where 
Vo =O tif, = 2.9150 + 7 1.3530 , (25a) 
when the function y = g, f, is used, and 
Yo =O +1 P, = 2.6914 + 7 1.3351 , (25b) 
when the function ® = G, F, is used. [Write in (24) G and H instead of g and h.] It is noted that? 
g(0) = 1, g(0)=0, (26a) 
h(O) == 0" 1) pe a (26b) 


The above eigenvalues (25a, b) are quite close to the correct one (13a), especially the one based 
on ®,. This encourages us to believe that the stress distribution derived from ®, will also be close 
to the correct stress distribution. The stresses derived from G, F’, were calculated and plotted in 1, 


and so were the end-face displacements 
Tr 


TAGs | [o, —» (09 +0,)]e-odr, (27a) 
w(t) = p> : [t0—» (6, +0,)].-0> (7b) 
wnt) = 75 { [o.—» (0, + 00)] dz, (27e) 
on) = som(0), Bis [ 2 ee eae a G11), (274) 


* The detailed expressions for the coefficients {, 8, and © are given in 1. For the sake of simpler notation, 
we shall frequently omit the subscript k of f, A, g, a, ete. 

* In general discussion, and in the specific case of normal boundary tractions, we use the symbols g andG 
for the z-function. In the specific case of loading by shear tractions we use the symbols h and H. Frequently 
Me pone eee! only about g or h when it is obvious, as in (26), that corresponding statements apply 
also to G an : 
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There is no direct way of judging the reliability of the stress results. But the displacements may 
be judged, and they were found to be disappointing. They were found to have the ranges 


(Oy = OY uy(1) = — 0.62, 
un(0) = 0, un(1) = — 0.03, 
w1(0) = 0.57, w1(1) = — 0.50, (28) 
; wrr(0) = 0.57, wy(1) = —0.01, 
wyn(0) = 0.57, wyu(1) = — 0.57. 


Equations (28) show that the various methods of displacement determination are quite inconsistent. 
One may regard the averages 

1 1 

= (ur + un) z (wy + wy + wit) (29) 


“ 


as perhaps representative of the true displacements; but the deviation of u, or u,,; from such an 
average is probably indicative also of the deviation of the calculated stresses from the true ones; 
and the indicated deviations are excessive, even in an engineering approximation. Instead of 
plotting similar graphs also for the function g,, or investigating, similarly, the higher modes, it 
looked more promising to preceed to the approximation B. 


5. Method B. We assume that 


le) = Ge). (30) 
For radial functions f,, F, we now choose 
f=kh, F, = (1 —k) F, (31) 
& ps 
LS 
29 BP 
LW 
28 
f = be hs 
0 
27 = 


08612 , 1.0 


Fig. 3. Variation, with factor k, of the eigenvalue «, + if, and the corresponding complementary 
energy U*, when the stress function pair p = k G,(z) f(r), ® = (1 — k) G,(z) F,(r) is used. 


where f,, F, are given by (21), (22). Then the boundary traction, see equations (6c, d), is still 

o,(r) = sr) =2°—1, (ih) 30 (32) 
when G(0) = 1, G’(0) = 0. In Fig. 3 of this paper are plotted the variation of U*, a5, B, with k. 
It is seen that the energy is minimum 


++” Ys, = 0.044626 ee) 


9 min 
at Resire 272 (33b) 
ia} 
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(Note also that the x, minimum occurs somewhat to the left of Uyj,.)1 One notices also, with 


considerable surprise, that the new eigenvalue 


Yo = 2.6564 + i 1.4067 (34) 


10 


—— 


-1.0 


Abb. 4a. 
Fie, 4a—b. End-face stresses and displacements, calculated by variational Method B, for loading of cylinder by end 
traction o°(r) = s,(r) = 2r?— 1. Uz, Uz], Wy, Wry, UzTT, Mo, Mp are based on Eqs. (27), (9). 


= 


Abb. 4b. 


1 A similar determination of the best k may also be made for shear loading b ul = (1/2 
ee lar d g by traction t, = (1/2) (r? —r). 
This determination has been carried out only for k = 0 and k = 1 d it 1 f a - 

= 0.005506, 0.006377, respectively. Ae EEE Se ae ae 
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is an inferior approximation to (13a) than is (25b), in spite of the fact the present 
solution is a much better solution — as judged by the compatibility of the displacements — 
than was the former one. In Fig. 4(a) we plot the end-face stresses as created by the traction 
distribution (32). We recall that the stress functions are now 


g=t G.f,. @® = (1— k)G, F,. (35) 


_where G is given by (24a), and k and y, by (33b), (34). In Fig. 4(b) we plot the displacements uy, 

UT, Wy, Wy, Wy. The great improvement over the corresponding curves of 1 for Method A is clearly 
indicated. We also show curves up, wy as obtained from the inapplicable formulas (9). [The latter 
curves should act as deterrents to the use of formulas of type (9) in problems where the functions 
are not strictly harmonic. ] 


6. Method C. We write 


(36) 


Then f, and F, are given by (21), (22), and the coefficients a;; etc., of Eqs. (18) are found to have 
the values listed below, in Eq. (37), in the first numerical column. The heading 10-? x indicates: 
multiply each entry by 10-*. The second column is based on the tractions s,, t, Eqs. (55). 


10-2 } 10-2 x 
Goo 35/3 27776/3 
Os 65/48 884/3 
ay, 7/64 5 
ay, —5/48 —68/3 
¢ 25/3 19840/3 
ae 65/48 884/3 (37) 
Coo 0,102284526 4.7207062 
Cos 0.236928104 59.1867160 
be) 85/3 27776/3 
bb, | (65/24 | 1768/3 
b,, |  0.158780410 | —4.9071620 
by, |  0.801334423 156.564392 
be | 3/8 136 


On taking the variation of U* one obtains 


Ep ee On Ok thas pr ue O8 Gee | ge 
ee OE) lie Og | 7g’ 88 I ag77 8 T OG I j 


bu d ox GSN x | of d On 5p RI 
a ~~ dz dg’ ' dz? ag!’ og J od dz + dg’ dz dg’’} © ae OH “2 0 
0 


Oe Apa? ; 
66-2 0G (38) 
+ |) ser 
from which there result the Euler equations 
2 [ag2 aM + (do2 — %1) vias + G98) + [bo. GY t+ (bo9 + boo by) CG" + by G] =9, (39a) 
[bo» én t (bog + doe b,3) 8 + ooog] + 2 [coo GIV + (¢y2 — €11) GC” + eG] =9. (39b) 
The eigenvalue equation is 
2 daa 94 + 2 (do2— M1) 7? + 2 Boo boo y* + (b20 4 bos by 1) y+ boo | = 6 (40) 
ean + (boo + Boe by1) 7? + boo 2 Cy ¥4 + 2 (Cog — O11)” + 2 “00 


and the solution is 
Yo = 2.664150 + t 1.401370 , (41 a) 


vv, = 4.919928 + i 0.5651619. (41b) 


V5 is a spurious root, resembling y,, whose purpose is to provide a better fit in the stress and dis- 
placement distributions than would be available on the basis of y, alone. 


lies 
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Note that by virtue of (6c, d) two boundary conditions have been implied in the statement 
of the problem of (36); namely, for the case of normal loading 

6, =S,{7),  t =": g(0) + G(0) =1, (42 a) 
g(0) +60) =0 (42b) 

and for the case of shear loading 
=0, 1 =t(r): (0) +H) =0, (43) 
il (43 b) 

These conditions, in turn, imply 
dg(0) = — 0G(0) , (44a) 
6g’ (0) = — 6G‘(0) . (44b) 


Introducing (44) into the last line of (38) we obtain the two supplementary (“natural”) boundary 


conditions 
og d a& epu d aX 


Meal Og’ dz og” 0G’ dz 0G’ ’ (45.3) 
AQ ag 
oS ee (45) 


Ze) 1 
In order to determine the coefficients in the axial functions 
go(z) = e—%* (m cos By z + nsin B, 2) + e—*2*(p cos f; 2 + q sin Bz) , (46) 
G,(z) = e~%2* (M cos B, z + N sin By z) + e—*2* (P cos B, z + Q sin fy 2) 


we substitute (46) into Eq. (39a). This provides four equations between the coefficients. Four more 
are obtained by substitution into (42a), (42b), (45a), and (45b). Solving the eight equations 


there result 


10 go(z) = e—%* [2.36108 cos Bz + 6.77826 sin fy z] + e—*2 [0.98573 cos fy z + 3.48640 sin fj 2] , 
(47a) 
10 G(z) = e—*2* [7.56523 cos By z + 1.25127 sin B, 2] + e—*22 [— 0.91204 cos B3 z— 3.88647 sin Bz) 
(A7b] 


for the normal loading case, and using (43a) and (43)), 
10 h,(z) = e—%* [0.94576 cos By z — 6.28375 sin B, 2] + e—92+ [2.35219 cos 6, 2 + 1.15598 sin fy z], 


(48 a) 
10 H,(z) = e—** [— 0.76619 cos 8,  — 12.92232 sin By z] + e—%9+ [— 2.53176 cos Bz — 1.38912 sin Bz] 
(48b) 
or the shear loading case. 
We may write the stresses in the form 
0, = [a,(r) cos B z + b,(r) sin B z] e—** + [e,(r) cos f’z + d,(r) sin f’ z] e—*'*, 
09 = [a9 cos + bo sin] e—** + [cg cos + dg sin] e—*'*, (49) 
c.=s() 8), G=g+6, 
tT = — t(r) G(z) 
One finds that for s,(r) loading 
10 a, = — 10.3718 r4 + 20.4374 r? — 10.0655, 
10 b, = 18.8499 r4 — 31.0229 7? + 12.1730 , 
10 c, = 0.5021 r? — 0.5161 r2 -- 0.0139 , (60a) 
10 d, = 1.1230 r4 — 3.1691 r2 + 2.0461 , 
10 a9 = — 5.3169 r4 + 14.7701 r2 — 10.0655 , 
10 bo = 7.6714 r4 — 6.4906 r? + 12.1730 , 
10 co = 0.5181 r4 + 0.4444 r2 + 0.0139 , (0B) 


10 dy = 10.1883 +4 — 14,0806 r? + 2.0461 
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and for f,(r) loading (S) = h + H taking place of ®) 


10 a, = — 5.9678 r* + 10.8119 72? — 4.8441, }- 
10 b, = 3.9372 r4 — 5.9586 r?2 + 2.0213 , 
10 c, = — 0.3179 r4 + 0.7336 r2 — 0.4157 , ea 
10 d, = — 0.1025 r4 + 0.4771 1? — 0.3686 , 
10 ag = — 2.7510 r*4 + 5.3479 r? — 4.8441 , | 
10 bp = 1.4327 r* + 0.3778 7? + 2.0213, 
10 co = — 2.1512 r4 + 2.7624 7? — 0.4157 ; kone) 
10 dg = — 1.7450 r4 + 2.6460 r? — 0.3686. 
The functions a(r),... are tabulated vs.r in Table 1, and the functions e—** cos f z, . .. are tabu- 


lated in Table 2. 


Table 1. The r-Factors of the Stresses (49), for s,, t, Tractions 


T | So | cr 5 d, a bo 77) | do 
| | | | 

0 —1.00 er —1.0066 | 1.2173 | 0.0014 | 0.2046 | —1.0066) 1.2173 0.0014 | 0.2046 
0.1 == al? a 9862 | 1.1865 © .0009 | .2014 | — .9918) 1.2108 0018 | 1906 
2 — .92 | — 9265 | 1.0962 | —-.0006:°| ..1921. | —~-.9483 | 1.1916 -0032 .1499 
3 Set Oe eek 9534 | — .0028 |}. .1770 | — .8779| 1.1602 0058 | 0861 
4 — .68 | — .7061 | 7692 | — .0056 |° .1568 | — .7838| 1.1177 0098 | — .0054 
5) == 5) |S | .0095 | — .0084 1324 | — .6705.| 1.0655 0157. | — .0837 
.6 — 28 | — 4052 3448. | — .0107 |; .1051 | — .5437] 1.0053 0241 | — .1702 
Bu! = 02 |) — 2541 1498. | — .0118 | .0763 » | —-.4105 9394 0356 | — .2407 
8 28 | — .1234 | — .0039 | — .0111 | .0478 “| — .2790 8704 .0510 | — .2792 
®) | 2 a= OSG 1.0588.) —— 0075.1, «10216. ° == 1590 .8012 .0714 — .2675 
899 -805 | — .0069 | — .0472 | — .0043 -0101 | — .1066 a MORE .0837 | — .2363 
1.0 OO 0 0) | 0 0 | — .0612 £7394 0976 | — .1846 


0 0 | —0.4844 0.2021 | —0.0416 | —0.0369 | —0.4844 | 0.2021 | —0.0416 |—0.0369 
0.1 —0.0495 | — .4734 .1962 | — .0408 | — .0364) — .4791| .2025 | — .0388 | — .0342 
2 — .0960 — .442] .1789 | — .0387 | — .0350 | — .4635| .2039 | — .0309 | — .0266 
53} P1309 ja 3919 1517 | — .0352. | — .0326 | — .4385|  .2067 | — .0184 | — .0145 
4 —= .1680')| — .3267 1169 = .0306 |— .0295 | = 4059"), > .2118—| — .0029- | .0010 
ae — L875) | — .2914)| 207 (Om =O 202020035 5019) ae. 2200) a 0140 | 0184 
6 — .1920 | — .1725 | 0386 | — .0193 | — .0210 | — .3275 |,  .2343 0300 | 0368 
ol = RES |) = AU .0047 | — .0132 | — .0160 | — .2884. 2550 | 04.21 .0509 
8 == 1440 |e — "0369 | — -.0180' |=—".0076 || — 0105 | —-.2948 .2850 0471 .0610 
ao) — .0855 | — .0002 | — .0222 | — .0030 | — .0049 | — .2317 .3267 0410 .0630 
£95 — .0463 .0053 | — .0149 | — .0013 | — .0022 | — .2258 3029 | 0325 0598 
1.0 Ona 0 0 0 | .0006 | — .2247 3002 | .0196 .0532 
Table 2. The z-Factors of th Stresses (49), for s,, t, Tractions 
z e “7 cos Bz ey OF sinBz ew 2 cos B' z oe *' % sin B'z (65) | G’ A) / 
| | | | 

0 1.0000 | 0 | 1.0000 | 0 } 1.0000 | 0 0 1.0000 
0.1 SS) * 0.1070 | 6104 | 0.0345 | 9626 | —0.6834 0.0764 .9067 
2 0040 | 1624 | 3714 0422 8741 | —1.0408 | mo 2072 
a) 4105 | Sige | 22093 0386 .1616 | —1.1802 | 1314 0632 
A .2918 a832au 1362 .0313 .6427 | —1.1809 |°  .1313 | — .0553 
Oe .1348 1507 | 0493 0174 A242 | — .9756 .1082 | — .1515 
8 -0516) 24] 1069 | 0176 ~ =. .0085 20do: 1 906943. | .0769 | — .1521 
1.0 0118 .0687 0062 | 0039 1440 | — .4473 .0495 | — .1195 
12 | — .0045 | 0406 | 0021 0017 0738 | — .2654 | .0294 | — .0824 
1.6 — .0088 AUTO 0002 | 0003 0126 |— .0725 | .0080 | — .0301 
2.0 — .0046 | 0016 | 0 0 — 10014 |— .0108 | .0012 | — .0078 
2.4 — .0016 — .0004 | 0 0 — .0023 .0023 | — .0003 | — .0010 
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The stress profiles o,(z, r), 0o(2, rT), 6,{r 2, ), T(z, r) are plotted, at z = 0, 0.1, 0.2, 0.4, 0.8, 1.2 vs. r, 
in Figs. 5(a) through 5(d) for the case of normal loading and in Figs. 6(a) through 6(d) for the case 
of shear loading. In Figs. 5(e) and 6(e) the variation, with z, of the stresses o,(2, 0) = oo(z, 0), 
oo(z, 1), o,(z, 0) = —o,(z, 1), and t(z, 0.6) is plotted. (¢ is a maximum at r = 1/)/3; at location 
r = 0.6, for which the calculations were made, t differs insignificantly from this maximum value.) 
The end-face displacements uy, uy, Wi, WH, Wr are plotted in Figs. 5(f) and 6(f).’ It is seen that, 
for sy loading, uy and uy as well as wy and wy coincide within graph paper accuracy. However, 
wy differs somewhat from wy and wy. Since (27d) involves — through du/dz — a stress derivative, 
the displacements wy, wy are vulnerable to errors in the stress derivatives, even if the stresses 
themselves are excellent approximations to the true stresses. Therefore, wy should be trusted in 
preference to wy, wi. 

Figures 5(f) and 6(f) show also the uy, wy curves, as derived from the inapplicable formulas (9). 
Comparison of Fig. 5(f) with Fig. 4(b) shows that even these curves are now greatly improved, but 
they still remain unsatisfactory. 


6. The rule of orthogonalization. In order to furnish approximations to the higher eigenvalues 
of Eq. (12) we must define an orthogonalization procedure for generating the higher polynomials 
Sa» F's; foo Pg; --- The approximate eigenvalue y, is obtained by placing f,, F), into the integrals (18) — 
and solving the corresponding secular equation (40). 


In 1 we were primarily concerned with the functions ®, = G, F,. These functions did not 
seem to be amenable to any “obvious” orthogonalization rule, so we defined the F;,(r) polynomials, 
of degree k +- 2, by the requirements that one coefficient in the polynomial be 1 (we have chosen | 
the r* coefficient for this purpose), two coefficients were determined by the requirements (la, lb) _ 
that t(z, 1) = o,(z, 1) = 0, and the remaining coefficients by the — for calculation purposes astro- 
nomically laborious — requirement that «, be an extremum. We called this method of construction 
a “generalized orthogonalization rule”’.2 ~ 
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Fig. 5a, Fig. 5b. 
Fig. 5a—f. Diagrams (a) through (e) illustrate the stresses Sys Oo Gyr T for loading of the cylinder 
by a°(r) = 8,(r) = 2r?—1. Diagram (f) compares the endface displacements u, w calculated by various methods [see Eqs. (7), (9)]. 


1 It is to be remembered that our u and w are denoted by 2 Gu, 2G w in conventional units. 


? See also, G. Horvay, Orthogonal Edge Polynomials in the Variational Soluti f 
Layer Problems in Elasticity, submitted to ZAMP. mmm ee otras 
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Fig. 5e. 


Tey (/) = + 0.011 
70,(/) = -040/ 
My (/) = -0.397 
My (/)= -0407 


Up (1)=-0.746 
U,(/) =-0.20/ 
Uy(/) = -0.203 


Fig. 5f. 


In the present case, where both the f,(r) and 
Fr) polynomials are used, the appropriate method 
of orthogonalization is, however, “obvious”. We 
found earlier that 


f,=>0—r), (52a) 
We now write 
fi(r) = 1 — 7°)? (ay + a le a 
few ee ay ig 1 *)s (52 b) 


This form of f,(r) guarantees satisfaction of 
the boundary conditions (la, 1b), and so these con- 
ditions need no longer be accounted for in an 
orthogonalization procedure. We therefore may 


define the orthonormalization rule for functions (52b) ina conventional fashion, by the requirement that 


[file fey dr = 


(53) 


‘ for k=1,)| 


0 fork+1.f 
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Fig. 6f. 
» 69, 9, t, for loading of the cylinder by 7°) = t,(r) = (1/2) @— r). 
Diagram (f) compares the endface displacements u, w calculated by various methods [see Eqs. (7), (9)]. 


One thus finds 


Fig. 6a-f. Diagrams (a) through (e) illustrate the stresses o, 


Si: = N, fi.» 
. 34 
f= = (Ie? , N= y10 : (54a) 


i ee 
fs = Gd—ryF(1—6r), N, = 40 yI4, 
(54) 
y 


it 
fs =gy(l— rR (1— 14? + 28 4), N, = 60 yI8. 
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Correspondingly, by formulas (5b, c), 
s, =~ (27 r!— 267° +4), 
= 


rs (225 1 — 315 rt + 114728), (55a) 


Ss = 


ty =-, Oro — 137 + 47), 


| 


(56 77 — 105.78 +. 577?—8r). (555) 


i 


ia 
on 


Boldface letters, as mentioned earlier, refer to the normalization (20). The normalization (53): 
(54a) then leads to 


sS; = k Sk » ek al Pa ip ON ins p= Ny Fe (56) 


It is noted with pleasant surprise that the characteristic tractions thus obtained differ but insignifi- 
cantly from the polynomials obtained in Fig. 3 of 1 by the method of extremizing «,. For in- 
stance, the requirement «, = max for Method A led in 1 to 


F,(r) = B, + B,r? + B,r* + B,r*; B,/B, = 0.34084 and 
Yq = 5.5888 + 1 2.4642 . 
In the present F', polynomial the coefficient ratio is B,/B, = 0.34450. The corresponding eigen- 
value calculated by “Method A” is 
Yq = 59.5832 + 7 2.4980. 


Now that we know f,,, s,, t,, we may obtain F, by integrating Eq. (5b). The explicit formula is 
1 u 


F,(r) = | faa a) | | oO) fi.(@) 7 Os — sty, | e? O—") fic(Q) de - (57) 


r 


In this fashion we find for instance that 


F,(r) = 2 [r> — 2.96914 r4 + 2.901961 r2 — 0.75127] . (58) 


The coefficients a;, etc. of (18) based on the functions f,, F, are listed in the second numerical 
column of (37). From (40) one obtains the eigenvalue 


74 = 5.433 + i 2.386 , 
yi, = 16.162 + 16.271. 


(59) 


Time was not available to calculate the stress and displacement distributions produced by trac- 
tions S,, {,. 


7. Stresses in a finite cylinder due to parabolic temperature distribution’. We return once more 
to the question of the stresses near the end z = 0 of an elastic cylinder (of length Z 4 radii) subject 
to the temperature distribution? 


O(z,r) =Oy,(1—r*). (60) 
The stresses far from the end-faces (indicated by subscripts A) are 
o4=—<(r—1), oo4=—(3—1), O4=52P—1), T4= 0 (61) 


1 In this section we revert to conventional stress units. ’ 
2 G. Horvay, Stress Reduction Obtainable in Sectioned Heat Generating Cylinders, Proc. Second Mid- 
western Conf. Solid Mech., Ann Arbor, Mich. (1956), p. 45. 
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On these there is to be superimposed the effect of end tractions 


6(r) =—= 2 =1). (62) 
The latter stresses! are — ©/2 times the stresses illustrated in Figs. 5(a) through 5(e). We now 
note a very surprising thing. The largest stress in the finite cylinder does not occur on the face 
z = 0. It is an axial stress [c. f., Fig. 5(e)] 
6(2.3, 1) = —9,(2.3, 0) = (05 eee 2) = 0.501 © (63) 
barely in excess of 6.4, max, at a relatively large distance z = 2.3 from the end-face. For the end- 
face z = 0 we obtain 


; wae 0.04. 
@.(0..0) = 4200.0) = e 0.25 + |e £3561) 059 © piaecg Oa ee (0.50 =e Al S=0486. 


| (64) 


The authors wish to express their appreciation to J. S. Gillespie, Manager of Engineering, 
Metallurgical Department of the General Electric Company, Detroit, Mich., for sponsoring, in part, 
this investigation, and for his interest. They wish to express their appreciation also to H. Poritsky 
and M. F. Bolster for permitting them to use their results, quoted in Eq. (13) and Fig. 2. These 
results are not available in the open literature. 


(Eingegangen am 29. August 1958.) 


Anschrift des Verfassers: Dr. G. Horvay, Schenectady, N. Y. (USA), Metallurgy and Ceramics Research 
Department, General Electric Company, General Electric Research Laboratory, Schenectady, N. Y. 


. 1 We denote these stresses by subscripts C; the resultant stresses are written without subscripts: 6 = o4 
+ oc. 

_ ® It will be recalled (p. 416 of Ref.) that in a thin-walled cylinder subject to a linear temperature gra- 
dient across the wall the largest stress is the end-face Og stress, 25 per cent higher than the on stress far from 
the end-faces. In Ref.® it was incorrectly estimated that also for the full cylinder the maximum stress is an 
end-face hoop stress, about 10 to 16 per cent higher than far from the end-faces, 

° S. Timoshenko and J. N. Goodier, Theory of Elasticity, New York 1950. 
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Large deflection of a circular cantilever beam 
with uniformly distributed load 


By K. Sato 


1. Introduction. Rohde has obtained the solutions for a straight cantilever beam with uniformly 
distributed load. Also Seames and Conway? have presented a special numerical procedure for 
calculating the large deflection of straight and curved beams subjected to forces concentrated or 
distributed. Although this method is simple and also gives comparatively exact results, it can 
be applied only to the case of non-inflectional elastica. 

The present paper deals with the bending of thin circular cantilever beam, convex downward, 
subjected to the uniformly distributed load. The basic assumptions are that the deformation is 
elastic and that the bending does not alter the length of beam. The deflections of the beam may 
be divided into two stages. In the first stage the deflections of the cantilever are relatively small 
and the curvatures at all points on the strained beam take the same signs as the initial ones, while 
in the second stage the strained beam has one point of inflection of the elastic line. 


2. The First Stage of the Deflections. Fig. 1 shows a circular cantilever beam under the action 
of a load of intensity w distributed uniformly along the arc length. It is assumed that the un- 
stressed form of the beam is a circular arc of radius R, as shown in dotted line in Fig. 1, and the 
included angle @ for arc is less than 77/2. 

Taking the origin at the free end B and 
denoting by @ the slope of the elastic line 
at any point of the strained beam, the 
Bernoulli-Euler equation may be expressed 
as 


dd M 1 
Piotart yy ice (1) 


where d@®/ds is the curvature at any point 
of the beam, M the corresponding bending 
moment, D the flexural rigidity, and ds is 
an element of the arc. Differentiating with 
respect to s and substituting (dM/ds) = 
(dM/dx) cos = w s cos M, we obtain 


iD Ws 
w= pe D. (2) 


Fig. 1. Illustration of the co-ordinates of a circular cantilever beam with 
uniformly distributed load, in the first stage. 


The end conditions in this stage are 


d@ 1 
at s=0: D=O,, de ae Re 
1 dP ®) 
a Reeth == Vs Sas SS, 
R s 
where I is the total length of the beam. We put 
Die) Ons. i? _ X'na,s'—!, is =e th ye Lae Ss , (4) 
n=0 ds nt ds? n=2 
where a, is the coefficient to be determined. The end conditions give aj = Pz and a, = — 1/R. 
Then equation (2) becomes 
se = > (cos P, cos T — sin Dz sin T), (5) 


F. V. Rohde, Q. Appl. Math. 11 (1953) p. 337. 
A eames and H. D. Conway, J. Appl. Mech. 24 (1957) p. 289. 


x 
WN 
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ie . . . . . 
where T' = 3/a, s". Expanding cos T and sin T and comparing the coefficients of various powers 


n=1 


of s on both sides, a, is determined as follows, including a) and a, obtained before: 


1 1 w 
G5 =D; Ca C= 05 a; = = 7, cos Dz , 
LS 2 at Ik My ay 
%=DR pin Oe: ds = —7 Fe pH 008 Oe > 
p/p \eae I a 
% = — 780 5) pie COO ean each 
IE PONG IE J au 5 lw DP 
PETER (>) (353 “jose ?;) + 7008 Rip 008 @s » 
maay AL of NB es SY 
a, — 6720 (5) sin Pz cos P, + 6720 R® D sin Pz 2 
w\s 1 Asia LW ae fie 19 Dir gets 
49 Eero) cos Dy [eieay ae0rn sin’ Dy) — 8 (D) (as5x0 — 799 Sin’ Op (6) 
| 1 1 w 
— $1840 Fp °° Pa» 
Ly SO) a thai lites 67 1 fw, 
a a(a) (ax3a00 sin Pz — 95599 Sin Dp) — 151200 Ri D) a 2 08. 
1 Loe. 
— 453600 RF DD Pe » 
1 /w\8 19 7 et ih pO Il es 9 
«= FD) cos Dp aan = SRg00-° | a (DB) (16800 739209 *12* Pp 
. 1 iG 
+ 7435200 7 Dp © a > 
Since dy/ds = sin ®, 
y = [sin @ds = T, sin Dy + T, cos Op , (7) 
é 
where ‘ 
FR po ~ "\ds = 1 3s ey ape as 58 No Eipiiae Se 
1 == |) ons 2, ons Sr Sa OS 5 41% — F59 P 
iy 
Gl ee 3 ~ OW ie =e 2 “ay (278) gh dae ae 
a peers) 2 as" Soairg Oy S° 7s | 5% 
Similarly, . 


* =f cosO@ds = T, cos Dp — T, sin By = >" gaara ot A : (8) 
9 n=3 (n — 2) = 


The bending moment M is obtained from equation (1) as 


d® 1 I ee 1 
M=D(Z+q)=D(Snas +3) (9) 


The slope at the free end of a beam is obtained from equation (4) and the end condition @ = 0 
atresia: 


oe oh a, 1", (10) 


In the limiting case R-> oo, these e 
straight beam obtained by Rhodel. 
flection of a circular cantilever bea 


quations for the first stage coincide with the results of the 


If the sign of w is reversed, these equations will give the de- 
m, concave downward. 


1 F, V. Rohde, Q. Appl. Math. 11 (1953) p. 337, 
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3. The Second Stage of the Deflections. Now, we proceed to the second stage of the deflections, 
in which the strained beam has a point of inflection. In this stage, it will be convenient to compute 
the deflection by reducing a circular cantilever to an initially straight cantilever with a couple M, 
applied at its free end. The couple M, must be constant, because it is the quantity relating to the 
initial shape of the beam, and can be expressed 


M,=—. (11) 


Then if we cut the strained beam at the point of inflection 0, as shown in Fig. 2, it may be treated 
as two straight cantilevers. The deflection of branch OA is equivalent to that of the straight can- 


Fig. 2. Illustration of the co-ordinates of a circular cantilever beam with uniformly distributed load, in the second stage. 


_tilever due to a distributed load w and a concentrated force ws, acting at the free end O, where sp 
is the length of branch OB of the beam. The deflection of branch OB is equivalent to that of the 
straight cantilever, of which the clamped end B is fixed with an angle ®, to the horizontal, sub- 
jected to a distributed load w and a concentrated force ws g acting at the free end O in the reversed 
direction to w, as shown in Fig. 2 (b). 

a) Branch OA. Taking the origin at the free end O and denoting the bending moment by M,, 
the Bernoulli-Euler equation d®/ds = M,/D gives 


dD 1 dM, w 
it Rees (s + sp) cos@. (12) 
Then putting 
ee d® oe 2D oe ; 
eae n oe eB n—|] IESE, ae NG! eae n—2 
Ce ne 5 7 he bes , 73 SA ob xStae (13) 


and denoting by ®, the slope at the free end O of the beam, the solutions to this branch are obtained 
in the same way as in the first stage. Coefficient b,, of power series is 


1 w 1 w 
2, , Decent. b, =—> pspcos PD, bs = — = 008 Ds | 
Bo —- SB B. . 
== tt 49, On ae )s 
where 
2 1 Nee 
i 0., B.=—3(5) Sp sin QD, cos D, , B=—4(5) sin D, cos Dy , 
w\3 1 Lae w\§ i ies 
= (5) sz cos Dy Cs — + sin? | , be (5) Sp cos Dy (T5 o- go sin?) : 


= 1 f(g w \4 ne es : 
_ ( B} cos D, (75 —s sint®,) + (5) sz cos Dy (a0 sin OD, — 50 sin? P| ; 
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Be =(5) cos cA et ae sin? P)— = 59 sim ‘| 
, w\! 4 
+ (3) Sp cos Dy ban sin D, — a sin’ Dp) : (14) 
By =— Ge) sp cos Dy ( aa -—— am sin? By + a sin'®,} 
“+ 3, cos Ds | e980 sin Dy a sin? Py) : 
dnse (5 ‘sh cos D l= sir sig 0500 sin? Dy sazi00 sin'®, 
= (3) sp cos Dy ( an sin D, =e GO sin? M, + — sin’ D 0) 
Sg = 3) Bee tase fie sna Sree: Again see 
mw (5) % re [a1 135200" sin Dy — at sin’ Py + 19950800 sin’ P, | 


Consequently the deflection Y, refering to X, Y co-ordinates (see Fig. 2), becomes 


Y = f sin @ ds = U, sin ®, + U, cos D , (15) 
where ‘ v 
i & i 1 1 
U, = [ cos( 3) 545") as aus = best 5 by bs 6 — (7,58 to baal as 
6 \ 
ea ae 1 1 ‘een ol I l 
ee [ sin (3b, s") a eh, pst tye a a by —7y ba)? Ming ae 
0 
Similarly - 
X = { cos Dds = U, cos ®, — U, sin D,. (16) 
f 


The slope of the elastic line at the free end O is 
Dy = — 3 by (Isp). (17) 


In this case the bending moment is 


However, it must be noticed that the bending moment M, is that of the straight cantilever beam. | 
Taking into consideration a couple My mentioned before and refering to x, y-coordinates, the bending | 
moment M in branch OA of the original circular cantilever is expressed as | 


1 co 
Me DR eC le 2) a , [2s 25,. (18) 


the bending moment in this case, we have 
@@ 1 dM’ 


ie Te (s — sp) cos@. (19) 
The end conditions in this branch are 
d®@ 
; Pati 23 ne } 
at s=0: Di= O,, Weert 
dp (20) 


| 
b) Branch OB. Taking the origin at the free end O as shown in Fig. 2(b) and denoting by M; 
As GS == Gar P= OT Back eee 
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For this branch, the solutions are obtained in an exactly similar way. Denoting by 6, the 
coefficient of power series for this branch, it is connected with the coefficient b,, for branch OA 
by the relations 


=—b, for odd n, 
(21) 


, 
b, —= 
, 
b= 65 for even n. 


The bending moment M, refered to x, y-coordinates, in branch OB of the original circular 
cantilever is 


M=D| >’ nb; (Gy —sy>* £4], r= sp. (22) 
n=2 
The end conditions at s = sp give Be 
Oz =D, + > bi sh, (23) 
n=2 
1 . ? on—l 
R= 2 bn Sh P (24) 


For convenience of the numerical calculation, we rewrite equation (24) in the form 


s TL es 1 : 7 i 
= Oo] a) 0s Py <r 7) V? cos Qpsin ® — V3 cos Dy ben — a sn 2) 
607 PAS) 
ey : i eT B40 0 y 
V4 cos ©, | Ta95510 sin Dy sey D,) + ; (25) 
3 3 
where Beer es 


By the proper connection of two branches at the point O, the whole shape of the strained beam 
are obtained. The procedure will be stated later. 

The shape of the strained beam having the point of inflection at its clamped end corresponds 

to that for the transition from the first stage to the second stage. Substituting sz = / and ®, = 0 

into equation (25), the relation between © and w /3/D in the transient state is obtained in the form 


l 2 , lwh fates Sl w [°\> 
© sey ee ee f— 2 | see 6 Go 
are PaaS 2D 2880 ( D } ™ 339500800 ( D ; GS) 
and from equation (23), denoting by ®z, the slope at the free end of a beam in this state, 
8 ay i ae Ill DN: 12787 w [3\> 
— [pee eer | a cliake areas y 
Px, Zab 3D * 22680 (" 4.670 265 600 ( D) a ey) 


4, Numerical Examples. The numerical calculations have been carried out for the cases 
O = 30°, 45°, and 60°. The deflection ratios (yy — yp)/l and (xg — x,)/l and the values of M1/D 
at the clamped end against the various values of wl?/D have been given for these three cases. 

It will be convenient to carry out the numerical calculations by the following procedure: For 
a given beam, the value of angle 0 = 1/R will be known. The values of wl*/D and ®,, in the 
transient state are easily obtained from equations (26) and (27), respectively, and the corresponding 
deflections can easily be determined using the equations for the second stage. 

; For the first stage, assuming a value of Dz (Og, < Og < O), a value of wl3/D is determined 
from equation (10). Using these values, equations (7) and (8) give the ratios y,/l and x,/I, respecti- 
vely. Then the deflection ratios (yy — yg)/l and (xg — %,)/l are determined. The non-dimensional 
quantity MI1/D at the clamped end is obtained from equation (9). 

For the second stage, assuming a value of ®, (0 << ®)) and solving equations (17) and (25) 
simultaneously, we obtain the ratio sp/l and the non-dimensional quantity wl?/D. Successively, 
equation (23) gives the value of ®g. Using these values, equations (15) and (16) give the ratios Y 4/I, 
Yz/l and X,/l, X%/l, respectively. Then the deflection ratios (y)— yg)/I [= (1 — cos @)/O + 
(Y4-+ Yp)/l] and (xg — %)/1 [= (X4 + Xz)/! — sin O/O] are determined. The non-dimensional 
quantity M1/D at the clamped end is obtained from equation (18). 

Repeating the calculations for various values of Dz, for the first stage, and ®), for the second 
stage, the deflection ratios and MI/D for the various values of wl?/D are found. These values were 
calculated with sufficient accuracy by using power series including eleven terms for the first stage 
and thirteen terms for the second stage. However, in the second stage also, there were many cases 


used eleven terms because good convergency of series. 
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Figs. 3 and 4 show plots of the deflection ratios and Ml/D at the clamped end versus wl*/D ; 
for 9 — 30°, 45°, and 60°. In Figs. 3 and 4, the marks o on the curves show the transition points 


from the first stage to the second stage. 
5. Summary. In this paper, the bending 


U of thin circular cantilever beam, convex down- 


Lé 


ward, under a uniformly distributed load is dis- 
cussedby the use of the Bernoulli-Euler equation. 
The solutions are obtained in the form of power 
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series. Numerical results are also presented. 
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Fig. 4. Variation of M1/D at the clamped end of a circular cantilever 
3 with wl?/D. 
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Fig. 3. Variation of (xg —x,)/l and (y)— yp)/! with wl4/D. 
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q raat Bezeichnungen und Parameter - Grundlagen der Schalentheorie: Geschichtliche Ubersicht. Biegetheorie 
der isotropen | Kreiszylinderschalen. Biegetheorie der Kreiszylinderschalen mit symmetrischen Ring- 
ras _ rippen. father der Kreiszylinderschalen mit Ringrippen unterhalb oder oberhalb der Schale - Flichen- 
asten. Membrantheorie : Geschichtliche Ubersicht. Das frei aufliegende Kreiszylinderrohr. Die Membran- 
theorie . . Isotrope Schalen mit Belastung am Lingsrand: Geschichtliche Ubersicht. Die charakteristische 
- Gleichung. Die Integrationskonstanten_ und die Multiplikatoren. Die Lésungsfunktionen und die Ver- 
= teilungsfunktionen fiir die Ringrichtung. Zahlentafeln fiir die Verteilungsfunktionen - Ringrippenschalen — 
mit Belastung am Langsrand : Geschichtliche Ubersicht. Schalen mit symmetrischen Ringrippen. Schalen 
mit exzentrischen Ringrippen - Schalen mit verinderlichem Halbmesser und veriinderlicher Schalenstiirke : 
-Geschichtliche Ubersicht. Vereinfachte Berechnungsgrundlage. Das Differenzverfahren - Berechnungs~ 
quicthoden fir Tonnendacher: Geschichtliche Ubersicht. Schalenberechnung nach der mathematischen 
Bra: Elastizitatstheorie. Allgemeine Balkentheorie, Modellschale. Tonnenreihe mit Randbalken ohne Biege- 
ee steifigkeit. -Tonnenreihe mit biegesteifem Randbalken. Einzeltonnen oder Aufenrand der Paralleltonnen. 
oe Der Einflu8 des Sekundarrandes und der Membrandeformationen, Streifenlasten P = sin Aé senkrecht 
zur Schalenfliiche . Trajektorienbewehrung, Binderberechnung : Geschichtliche Ubersicht. Die Trajektorien- 
bewehrung. Binderberechnung - Stabilitit, Ring- und Lingsmomente: Geschichtliche Uhersicht. Eigen- 
werte des. Ringknickens. Eigenwerte des Axialknickens. Bemessung der Schalendruckzone, Ring- und 
= Langsmomente - Das vorgespannte Tonnendach: Geschichtliche Ubersicht. Berechnungsgrundlagen. Vor- 
: : gespannte J rremcearaes _ Vorgespannte Einzeltonne mit frei aufliegenden Randbalken » Belastung am _ 
Ringrad x0; Geschichtliche Ubersicht. Die Lésungsfunktion. SchnittgréBen und _Verschiebungen. 
_Voragautine Literaturverzeichnis - Sle a ant 


: i 


is Die Berechnung der weitgespannten Tylicdecechalon kann ech der ieaestied Balkentheorie erfolgen. 
Bei mittleren und Ste Spannweiten versagt dieses Verfahren, und man ist auf elastizitatstheoretische 
Berechnungen angewiesen. Zahlentafeln erleichtern die notwendigen Berechnungen. Trotz der Verein- 
-fachungen, die dadurch erzielt werden kénnen, bleibt die Berechnung der kleinen und einfachen Schalen 
noch zu zeitraubend. Eine praktisch brauchbare Lisung ist, die Ringverteilung nach der Elastizitats- 
_ theorie und die Langsverteilung nach der elementaren Balkentheorie zu berechnen. Dieses Buch geht 
einige Schritte weiter auf diesem Wege. ; Es wird eine ,,Modellschale“ eingefiihrt, und diese Schale wird 
so gewahlt, daB erstens genaues Berechnen und Tabulieren der Schnittkrafte méglich ist, und zweitens 
80, daB die Randbedingungen bestméglichst mit der gegebenen Schale tibereinstimmen. Es werden be- 
stimmte ‘Verhiltnisse zwischen den Randkraften gewahlt, und die Randglieder werden den gewahlten 
Verhiltnissen angepaBt. Zahlentafeln einiger Modellschalen mit und ohne Ringrippen werden gegeben. 
Mittels dieser Zahlentafeln 1aBt sich die Spannungsverteilung schneller berechnen als nach der elementaren 
Balkentheorie, auch bei groBen Spannweiten. MaBgebend fiir die Betonabmessungen ist die Stabilitat 
der Schale. Die Knickberechnungen sind darum eingehend untersucht, und ein allgemeines Verfahren 
ist ausgearbeitet, das sowohl fiir Stabe, Bogen, Platten als auch fiir Schalen gilt. Die genaue Grundlage 
: der Schalentheorie wird untersucht und es wird gezeigt, wie eine iterative Berechnung einfach und mit 
beliebiger Genauigkeit durchgefiihrt werden kann, Zahlentafeln der Integrationskonstanten und Multi- 


plikatoren sind beigegeben. 
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INHALTSUBERSICHT: 


Einiges iiber die Grundlagen der Baustatik: Einfiihrung. Zusammensetzung und Zerlegung von Kraften. 
Darstellung von Auflagerungsméglichkeiten. Bestimmung des Grades der statischen Unbestimmtheit. 


Berechnung von Trigheitsmoment, Zentrifugalmoment und Schwerpunkt. Parabelkonstruktion. 


Das Verfahren der Kinematik: Allgemein. Die Bewegung der starren Scheibe. Die zwanglaufige kine- 
matische Kette. Polplane. EinfluBlinien mit Hilfe von Polplanen. Feststellung der Verschieblichkeit 
mit Hilfe von Polplanen. Das Verfahren der F'-Figuren. Feststellung der Verschieblichkeit mit Hilfe der _ 
F’-Figur. Berechnung paver Gchuittl watt ais Tite ieee -Figur, c 


Formveranderungen: Die elastischen Formanderungen. Formanderungsarbeit. Das Prinzip der virtuellen 
Verriickungen und der virtuellen Arbeiten bei elastischen Kérpern. Die Satze von der Gegenseitigkeit 
der elastischen Formanderungen, Satz von CASTIGLIANO. Die Mohrschen Satze. Auswertung von Form- 
anderungsintegralen. Das Verfahren der w-Zahlen. Verfahren der Wm-Gewichte. Der Williotsche Ver- 


schiebungsplan. Formanderungen an statisch bestimmten Tragwerken. 


Verfahren zur Auflésung linearer Gleichungen: Auflésung mit Hilfe von Determinanten. Eliminations- 
verfahren. ; 


Berechnung statisch unbestimmter Systeme: Kraftgré8enverfahren. Formanderungsverfahren. Iterations- 
verfahren. Kombinierte Verfahren. 


Sonderkapitel der Statik - Aufgaben - Hilfstafeln - Sachverzeichnis. 


Sinn und Zweck dieses umfassenden Buches tiber Baustatik ist, dem Studierenden und dem in der Praxis 


tatigen Ingenieur zu helfen. 


Die Theorien der wichtigsten Verfahren werden in angemessenem Umfang in Jeichtverstandlicher Form 

behandelt und ihre Anwendung jeweils durch eine Anzahl sorgfaltig ausgewahlter, vollstandig durch- 2 
gerechneter Zahlenbeispiele gezeigt. Oftmals sind dieselben Beispiele nach verschiedenen Methoden 
berechnet worden, so daf der Leser ein Gefiihl fiir die ZweckmaBigkeit des zu wahlenden Verfahrens 
bekommt. = 
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